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THERMOSPHERIC  PROPAGATION  OF  SONIC  BOOMS 
FROM  THE  CONCORDE  SUPERSONIC  TRANSPORT 

Executive  Summary 

In  recent  months  accounts  have  arisen  which  attribute  at  least  some  of  the  so  called  "Mys- 
terious East  Coast  Booms"  to  operation  of  the  Concorde  aircraft.1  2 These  reports  differ  from 
the  conclusion  of  the  Naval  Research  Laboratory  investigation3  which  found  that,  while  the 
Concorde  was  probably  responsible  for  the  booms  heard  in  Nova  Scotia,  the  East  Coast  booms 
were  most  likely  caused  by  high  performance  military  aircraft  operating  in  the  warning  areas  off 
the  East  Coast.  Conjectures  about  Concorde  induced  booms  have  been  encouraged  by  theories 
advanced  by  Dr.  Jeremy  Stone  of  FAS  and  Dr.  Richard  Garwin  of  IBM4  involving  long  range 
propagation  of  the  Concorde  sonic  boom  through  the  thermosphere,  production  of  superbooms, 
and  possible  adverse  effects  on  the  thermosphere.  These  conjectures  have  raised  some  serious 
questions  about  the  nature  of  sonic  boom  propagation  in  the  upper  atmosphere.  That  a ther- 
mospheric mode  of  propagation  exists  is  not  in  question.  Evidence  for  such  a propagation 
mode  already  exists  and  is  well  known.5  However,  the  speculations  were  based  on  an  over 
simplified  model  of  sonic  boom  propagation  and  lead  to  a number  of  erroneous  conclusions. 

In  response  to  a request  by  Dr.  Frank  Press,  Director  of  the  Office  of  Science  and  Tech- 
nology Policy,  NRL  has  continued  its  investigation  and  in  particular  looked  in  detail  at  Dr. 
Garwin’s  suggestions.  In  this  report  we  describe  a model  which  allows  an  accurate  description 
of  the  effect  produced  by  the  portion  of  the  Concorde  sonic  boom  which  reaches  the  thermo- 
sphere. We  reach  a number  of  conclusions  which  differ  considerably  from  those  proposed  by 
Dr.  Garwin. 

Scope  Of  The  Investigation 

We  consider  only  those  sonic  booms  which  are  refracted  to  the  ground  by  the  sound  velo- 
city gradient  in  the  thermosphere.  That  is,  we  consider  only  those  waves  which  reach  an  alti- 
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tude  of  more  than  100  km  before  curving  downward  and  returning  to  earth.  We  consider  both 
the  sonic  boom  which  propagates  directly  to  the  thermosphere  and  the  boom  which,  while  ini- 
tially traveling  downward,  reflects  off  the  water  and  returns  to  the  thermosphere.  We  recognize 
that  off  the  flight  track  the  sonic  booms  can  also  be  returned  to  the  ground  by  winds  in  the  stra- 
tosphere, usually  from  a height  of  about  50  km.  These  stratospheric  returns  have  been  studied 
by  the  British6  and  are  believed  to  be  the  cause  of  the  "bumps  in  the  night"  heard  in  England. 
This  propagation  mode  is  well  understood  and  shall  not  be  considered  in  this  report. 

We  use  a realistic  model  for  the  atmosphere  but  do  not  take  into  account  the  details  of 
the  wind-temperature  structure.  The  reason  for  limiting  the  problem  in  this  manner  is  that  the 
atmospheric  structure  depends  on  the  time  of  day,  the  season,  the  state  of  solar  activity,  lati- 
tude, weather  etc.,7  the  variability  of  which  render  consideration  of  the  detailed  structure  practi- 
cally meaningless.  We  rather  aim  towards  obtaining  both  typical  and  worst  case  results. 

We  consider  only  the  case  of  the  aircraft  in  straight,  level,  and  steady  flight.  We  recog- 
nize that  it  may  be  possible  for  a upward  going  sonic  boom  to  focus  at  the  ground  for  a 
decelerating  or  turning  aircraft.  Such  a focus,  however,  would  occur  over  only  a small  area  on 
the  ground  and,  in  any  case,  would  result  in  ground  pressure  levels  which  would  exceed  our 
predictions  by  no  more  than  a factor  of  three. 

We  include  in  our  theory  the  variation  with  altitude  of  sound  speed,  density,  and  pressure 

resulting  from  a realistic  model  of  the  atmosphere.  We  also  include  the  nonlinear  stretching 

/ 

and  decay  of  the  wave  the  focussing  which  occurs  at  the  caustic  near  the  turning  point,  the 
effects  of  linear  acoustic  attenuation,  and  the  effects  of  winds.  From  our  model,  we  determine 
the  boom  strength  as  a function  of  altitude  and  the  ground  pressure  both  on  and  off  the  flight 
track.  We  have  studied,  and  determined  to  be  unimportant,  gravity  wave  effects,  strong  shock 
effects,  plasma  effects,  and  the  effects  of  a finite  mean  free  path. 


Results  And  Conclusions 

We  have  developed  a model  for  thermospheric  propagation  of  the  sonic  boom  from  a 
Concorde  aircraft.  The  model,  which  is  for  the  most  part  analytic,  accounts  for  all  the  physical 
effects  which  we  believe  to  produce  a significant  effect  on  the  result.  Based  on  the  results 

% 

obtained  from  this  model,  we  reached  the  following  conclusions: 

1.  The  ratio  of  boom  pressure  to  undisturbed  pressure  in  the  thermosphere  never  exceeds 
0.2.  This  confirms  the  validity  of  our  weak  shock  model  and  implies  that  no  large  perturbations 
(such  as  high  temperatures  which  might  produce  strange  lights)  of  the  thermosphere  will  occur. 

2.  Over  90%  of  the  sonic  boom’s  energy  is  lost  before  it  reaches  100  km.  There  is,  there- 
fore, insufficient  energy  remaining  in  the  wave  to  cause  winds  nearly  as  large  as  those  which 
naturally  occur  in  the  thermosphere. 

3.  The  sonic  boom  will  reach  the  ground  at  a horizontal  distance  of  about  320  km  from 
the  point  it  left  the  aircraft. 

4.  Along  the  flight  track,  the  Concorde,  if  it  remains  in  straight  and  level  flight,  will  pass 
over  head  approximately  13  minutes  before  the  first  arrival  of  the  boom  from  the  thermo- 
sphere. The  average  horizontal  velocity  of  the  sonic  boom  is  less  than  half  that  of  the  Con- 
corde in  cruise  mode.  It  is  thus  highly  unlikely  that  there  could  exist  reasonable  conditions  of 
aircraft  deceleration  under  which  the  boom  would  precede  the  aircraft. 

5.  The  sound  received  along  the  flight  track  is  a minimum.  The  largest  signals  occur  200 
to  400  km  from  the  flight  track. 

6.  The  largest  signals  are  produced  by  the  initially  downward  boom  which  reaches  the 
thermosphere  by  first  reflecting  from  the  ocean. 
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7.  The  sound  received  will  in  all  cases  be  small  (less  than  .5  Pa)  and  very  low  in  fre- 
quency (approximately  .1  Hz). 


We  conclude  that  thermospheric  returns  from  Concorde  are  of  sufficiently  low  amplitude 
and  frequency  that  it  is  unlikely  that  they  are  either  responsible  for  the  East  Coast  events  or 
likely  to  disturb  the  public.  In  addition,  we  conclude  that  the  sonic  boom  from  the  Concorde 
does  not  have  sufficient  amplitude  or  energy  to  produce  a deleterious  effect  on  the  thermo- 
sphere. 


1.  Geometric  Acoustics,  Snell’s  Law  and 


Blokhintsev  Invariance 

The  propagation  of  weak  sound  waves  in  the  atmosphere  under  a limited  number  of 
assumptions  is  governed  by  the  theory  of  geometric  acoustics.  This  is  a special  case  of  the 
more  general  geometric  theory  of  general  linear  wave  propagation.  A general  theory  for  non- 
dissipative  waves  may  be  derived  (see  for  instance  Hayes8)  which  is  an  asymptotic  theory  for  a 
slowly  varying,  nonuniform  medium  for  which  the  local  solutions  are  very  close  to  those 
obtained  for  plane  waves.  The  theory  is  asymptotic  in  the  sense  that  it  recaires  the  wave 
number  to  be  large  in  some  relative  sense.  For  example,  the  wave  length  must  be  small  com- 
pared to  the  scale  length  of  the  inhomogeneities  in  the  atmosphere.  If,  in  addition,  the  waves 
are  nondispersive,  arbitrary  functions  of  some  phase  variable  may  be  used  to  describe  a general 
wave  signature.  A linear,  inviscid,  acoustic  theory  leads  to  a geometric  theory  which  is  nondis- 
sipative  and  nondispersive.  The  geometric  theory  allows  the  partial  differential  equation  for  the 
wave  fronts  to  be  reduced  to  the  integration  of  a set  of  ordinary  differential  equations  along  the 
characteristics.  It  also  provides  a means  of  describing  the  kinematics  of  the  solution  indepen- 
dent of  the  dynamics.  Thus,  from  different  orders  of  the  asymptotic  expansion,  we  obtain  solu- 
tions for  the  ray  paths  and  wave  front  curvature  and  an  energy  invariant  rule  which  allows  the 
calculation  of  the  wave  amplitude.  The  equations  for  the  ray  path  and  wave  front  vector  given 
for  geometric  acoustics  are:  (see  for  example  Hayes  et  al).9 


— « c ” cn  + u (1-1) 

dt 

— Vc  + (Vu)n  — nlnVc  + n(Vt7)n]  (1-2) 

dt 


where  7 is  the  ray,  c the  local  sound  speed,  u the  wind  vector,  and  n the  wave  front  normal. 
For  the  case  of  a stratified  medium  this  equation  possesses  a first  integral  which  yields  Snell’s 
law  (see  Ref.  9).  In  simple  terms,  this  law  states  that  the  horizontal  vector  component  of  the 

inverse  phase  velocity  — is  a constant  along  rays.  That  is. 


where 


cos  9 n 
Cn 


const. 


n = cos  tin'  4-  sin  tik 
c„  = c-ii  = c + u-n 


(1-3) 


The  second  order  term  of  the  expansion  in  reciprocal  wave  number  leads  to  an  energy 
invariance  law.  This  is  equivalent  to  the  definition  of  an  adiabatic  invariant  for  the  wave  propa- 
gation. Restricted  to  the  validity  of  geometric  acoustics,  an  invariant  law  was  found  by  Blo- 
khintsev,  which  is  valid  for  an  atmosphere  with  steady  winds10.  This  result  has  been  extended 
to  be  valid  for  unsteady  medium  by  Garret"  and  Hayes,12  but  this  has  little  application  to  the 
problem  at  hand  and  is  difficult  to  implement.  The  basic  result  is  that  an  invariant  is  found 
which  is  the  product  of  the  Rayleigh  acoustic  energy  density  pq 2 and  the  ratio  of  the  local  phase 
velocity  to  the  local  sound  speed  c„/c.  The  product  of  the  invariant  quantity  times  the  volume 
flow  is  then  a constant  along  a kinematic  ray  defined  earlier. 


C — 

pq2  — c ■ A - const.  (1-4) 

c 

We  choose  to  express  this  invariance  following  Hayes  et  al4  as  a constant  c'  times  the  square  of 
a variable  V 


c ' y2  — ~ cnc-/4  (1-5) 

c 

c is  the  group  velocity  defined  earlier  and  A is  any  cross  section  which  completely  cuts  the  ray 
lube.  This  then  provides  the  basis  by  which  the  amplitude  of  the  wave  can  be  computed  along 
a ray  path  knowing  the  amplitude  anywhere  along  the  ray  path.  If,  as  is  convenient  for  a hor- 
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izontally  stratified  atmosphere,  we  take  A to  be  the  area  cut  by  a plane  of  constant  r,  then 
c-.4  = c-  4.  From  Snell's  law  the  result  can  be  obtained  that  c.  - c sin  9 and  c„  - c'  cos  8. 
This  gives  our  conservation  variable  as: 


(1-6) 


' I = y/ftt''A  sin  8 cos  8 

The  disturbance  Mach  number  (q/c)  is  a measure  of  the  strength  of  the  signal.  The  problem  is 
completely  determined  once  the  ray  tube  area  A is  computed.  This  can  be  done  by  quadrature 
once  the  solution  for  the  rays  is  known.”  12  The  theory  as  it  stands  is  a completely  linear  one. 
This  leads  to  physically  incorrect  results  as  we  shall  see  later. 


2.  Linearized  Lossless  Approach  (Garwin)4 


As  long  as  the  linear  theory  of  geometric  acoustics  remains  valid,  the  Blokhintsev  invari- 
ance allows  signal  strength  to  be  calculated  for  an  arbitrarily  varying  inhomogeneous  atmo- 
sphere as  long  as  it  remains  steady.  We  may  use  this  invariance  to  estimate  the  signal  strength 
for  a supersonic  aircraft  such  as  the  Concorde  at  the  altitude  where  the  rays  turn  over.  For 
simplicity,  let  us  assume  there  are  no  winds  and  assume  a simple  cylindrical  spreading  of  the 
wave  so  that  the  ray  tube  area  is  proportional  to  the  distance  travelled.  It  can  be  shown  for  a 
horizontally  stratified  atmosphere  that  this  distance  is  the  horizontal  distance  if  we  are  looking 
along  the  flight  track.  We  are  neglecting  the  focussing  of  the  wave  near  the  caustic,  an 
unnecessary  complication  at  this  stage.  The  invariant  is  then 


c’  V 2 = pq 2 cA„  = — A„  = const 
pc 


(2-1) 


where  bp  = pcq  when  geometric  acoustics  is  valid.  If  we  use  the  atmosphere  described  in  the 
model  atmosphere  section,  we  find  that  for  an  M — 2 aircraft  at  17.5  km,  the  altitude  at  which 


the  phase  velocity  equals  the  aircraft  speed  is  about  160  km.  The  ray  path  found  from  using 
the  model  atmosphere  extends  about  150  km  horizontally  to  the  turnover  point.  If  we  use  Dr. 
Garwin’s  estimates  for  the  boom  strength  of  100  Pa  at  25  km,  we  can  calculate  the  boom 
strength  at  160  km.  However,  the  more  interesting  number  is  the  ratio  of  the  acoustic  energy 
to  the  local  ambient  energy  density  since  this  is  a measure  of  relative  boom  strength.  To  get 
this,  we  observe  that 


Thus,  the  ratio  of  (bp/p)2 


A,.  **  pcA„  = const. 
pc  yp2 

at  160  km  compared  to  the  value  on  the  ground  is  given 


by 


(2-2) 
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(A p/p)  jto  ^ (pcA ) y ^ | P«  | |_£*_j  | (2-3) 

(A p/p)\  (.PC A)  iso  I P 160  ) | c 160  J | A 160 

The  ratios  for  this  example  are  given  by:  pjp\w  ” 3.7  x 108,  and  cK/c\w " 0.5,  where 

,4„/.4  IM  - ~ - 1/6,  yielding  a ratio  of  - 3 x 107.  This  represents  a very  large 

* 150  (A p/p); 

increase  in  the  strength  of  the  wave  which  would  completely  invalidate  the  geometrical  acoustic 

theory.  As  we  shall  see  in  the  next  section,  there  are  other  factors  which  become  important 

rendering  the  above  simple  estimate  invalid. 


This  result  is  similar  to  Dr.  Garwin’s  argument4  when  he  postulated  strong  shocks  and 
large  temperature  rises  in  the  thermosphere.  In  his  approximation  he  used 

AJA  uo  - - 05  but  assumed  cv/c,80  - 1 Also  he  used  pjp ,80  - 109  giving  a result  that 

(\p/p)  iso  _ 5 x jq7  wj,jch  differs  little  from  the  estimate  given  above.  The  critical  assump- 
(A p/p)  l 

tion  in  the  above  estimates  is  that  the  wave  will  propagate  to  the  thermosphere  as  a linear  wave. 
This  is  not  true  even  for  weak  waves  as  we  shall  see  in  the  next  section. 


3.  Nonlinear  Analysis 


The  analysis  in  the  previous  section  was  based  on  a linear  geometrical  acoustic  theory 
used  to  derive  ray  equations  and  an  energy  invariant  law.  It  is  well  known,  however,  that 
acoustic  signals  distort  when  propagating  over  long  distances.  This  distortion  is  caused  by  a 
weak,  nonlinear  effect  resulting  from  small  changes  in  propagation  speed  which  are  proportional 
to  the  local  strength  of  the  disturbance.  Although  the  nonlinear  effect  is  locally  weak  (an 
assumption  that  must  be  true  for  geometric  acoustics  to  be  valid),  it  is  cumulative  and  the  dis- 
tortion of  a signal  which  has  propagated  over  a long  distance  may  be  far  from  small. 


In  order  to  simplify  the  analysis,  we  make  the  usual  assumption  that  to  lowest  order  the 
phase  shift  of  the  wave  due  to  the  change  in  propagation  speed  is  the  only  nonlinear  effect  con- 
sidered. Hayes13  has  shown  that  for  the  case  of  level  flight  in  a uniform  atmosphere,  the  error 
corresponds  to  a higher  order  term  in  a small  parameter.  The  assumption  of  a single  nonlinear 
effect  should  be  valid  as  long  as  geometric  acoustics  is  valid  and  as  long  as  the  disturbances  are 
weak.  The  problem  of  the  nonlinear  effect  on  rays  has  been  considered  by  Whitham14  15  16  for 
finite  strength  shocks. 


The  effect  of  the  small  disturbance  on  the  propagation  velocity  is  to  change  the  group 
velocity  from  c - cn  + /7  in  the  undisturbed  fluid  to  c = (c  + q)n  4-  U.  The  change  in  the 
speed  of  sound  is  given  by  Ac  = (T  — 1)  = (T  — 1)<7  where  I is  a thermodynamic  quantity 


defined  by  F - 1 = 
shift  expressed  as  a 


c 


etc 


. The  change  in  propagation  velocity  is  given  by  Yqn. 


distance  from  the  linear  wave  front  is  given  by 


The  phase 
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r 


d!  r 

J<--rQ 

where  l is  the  time  along  the  ray  and  / is  the  local  displacement.  The  local  displacement  is 
related  to  the  phase  by  dl  — c„dr  where  r is  a phase  variable.  Thus. 

(3-1) 

dt  c„ 

Integrating  Eq.  (3-1)  yields  the  phase  shift  of  the  nonlinear  signal. 


The  linear  theory  gives  the  result  that  for  the  Blokhintsev  invariant 


That  is,  along  the  linear  characteristics,  the  quantity  V — constant.  If  we  define  a phase  variable 

dr  - dt  — — where  din  - dr,  then  the  effect  of  the  nonlinearity  can  be  included  with  an  extra 

term  proportional  to  the  phase  shift  giving  an  inviscid  Burgers’  equation. 

^ = la  ^ (3-3) 

dl  c„  dr 


If  we  now  use  the  definition  of  V to  relate  it  to  — and  if  we  introduce  an  "age"  variable 

c„ 

defined  by 


c„(pA  sin  8 cos  9)U1 

then  the  equation  becomes  a simple  expression  with  the  independent  variables  "age"  and 


retarded  time  (or  phase).  Thus, 


d£  _ v 8K 


df  dr 

This  equation  has  an  analytic  solution  in  terms  of  the  linear  phase  variable  t,.  The  phase  shift 

is  given  by  ~~  - - K(r,  O-  Now  if  we  assume  that  the  Blokhintsev  invariance  and  ray  tube 
df 

area  are  unaffected  by  the  phase  shift  (an  assumption  we  justified  earlier),  then  we  have 
V(t,  f)  — K (r |)  where  t*  is  the  linear  phase  and  r the  actual  phase.  The  solution  of  the 
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phase  shift  equation  gives  r = rj  - £ F(t,).  The  linear  phase  r,  is  then  a function  of  t and 
The  solution  of  this  equation  may  become  multivalued  in  t.  This  can  be  resolved  by  taking 
into  account  the  presence  of  shock  waves. 

The  distortion  of  an  arbitrarily  shaped  signal  is  depicted  in  Figure  3-1.  The  rate  of  distor- 
tion is  proportional  to  the  amplitude  of  the  wave  for  any  value  of  the  phase  variable.  Assuming 
Blokhintsev  invariance,  the  amplitude  of  the  signal  remains  unchanged.  Thus,  the  portion  of 
the  signal  with  large  amplitude  will  tend  to  overrun  the  region  ahead  of  it  with  smaller  ampli- 
tude until  the  solution  becomes  multivalued.  This  multivaluedness  may  be  removed  by  fitting 
shocks  to  the  solution  rendering  it  single  valued.  For  weak  shocks  it  can  be  shown  that  the 
shock  speed  is  equal  to  one  half  the  sum  of  the  characteristic  speed  in  front  of  and  behind  the 
shock.17  This  leads  to  the  relation 

— ( V.  + VJ  <3-6) 

^ SHOCK  ^ 

This  can  also  be  shown  to  be  equivalent  to  the  equal  area  rule18  which  says  that  the  integral  of 
the  solution  with  respect  to  phase  r with  the  shock  fitted  should  give  the  same  area  as  the  mul- 
tivalued solution  for  each  lobe  of  the  solution.  This  is  also  equivalent  to  assuming  the 
existence  of  a velocity  potential  which  is  continuous  across  shocks.  Again  shocks  are  located  so 
that  the  integral  of  the  wave  with  and  without  shocks  gives  the  same  result. 

As  shown  in  Figure  3-1,  as  the  signal  distorts  further,  shocks  tend  to  overrun  one  another 
and  eat  up  the  entire  solution  until  the  asymptotic  form  of  the  pulse  forms  an  A wave  (which  is 
characteristic  of  the  far  field  solution)  with  two  shocks  separated  by  a linear  expansion  region. 
This  form  of  the  solution  will  always  be  arrived  at  for  a sufficiently  large  "age"  parameter. 

In  the  case  of  an  A wave,  there  exists  a particularly  simple  analytic  solution  to  the  fitting 
of  shock  waves.  The  solution  for  the  linear  part  of  the  wave  is  given  by 


1 + Hi 
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The  bisector  rule  gives  the  velocity  of  the  leading  shock  as: 


dr 

dZ 


I SHOCK 


This  can  be  integrated  to  give 


1 y - l 

2 5 2 1 + fit 


(3-8) 


t«  = r „(1  +M1'2 

From  the  initial  conditions  at  £ = 0,  we  get 

This  solution  gives,  as  required,  that  the  area  under  the  front  half  lobe  is  conserved. 

Kr„(l  + «) 


2 J Vdr  = F,r,  = 


= V T 

1 + 


(3-9) 


(3-10) 


(3-11) 


The  error  involved  in  using  the  far  field  N vr  c solution  for  the  entire  path  involves  only 
the  period  of  aging  the  wave  into  an  N profile.  For  cases  of  upward  propagation,  where  the  age 
becomes  very  large,  this  takes  place  rather  quickly.  The  error  involved  would  be  a somewhat 
enhanced  dissipation  of  the  wave  by  shocks  in  the  early  part  of  the  calculation.  An  equivalent 
N wave  for  a given  F function  for  the  aircraft  can  always  be  obtained. 


In  order  to  complete  the  solution  to  the  inviscid  Burgers’  equation,  we  must  evaluate  the 
"age"  integral  for  £ as  a function  of  position  along  the  ray  paths.  Recall  that 

Yd! 


c„  (pA  sin  9 cos  9) 


1/2 


(3-12) 


From  Snell’s  law  we  have  that 


dz  „ 

— = c sin  9 
di 

c„  = c‘  cos  9 


(3-13) 


Initially,  so  that  the  integrals  will  be  reasonably  simple  and  analytically  tractable,  we  make 
a few  simplifying  assumptions.  One  is  that  there  are  no  winds,  and  the  second  is  that  only  rays 


along  the  flight  track  will  be  considered.  Later  we  will  dicsuss  how  to  remove  these  restrictions 
at  the  expense  of  complicating  the  integrals.  With  no  winds,  c = c'  cos  0 and 


r r dz 

c'2  ^ sin  9 cos2  9 (pA  sin  9 cos  9)U2 


(3-14) 


Our  model  atmosphere,  discussed  in  section  four,  has  postulated  two  separate  parts  to  the  ray 
path.  In  the  lower  region,  up  to  a break  point  altitude,  the  sound  speed  is  held  constant  yield- 
ing straight  ray  paths  at  the  initial  inclination  angle.  The  second  part  consists  of  a region  where 
the  sound  speed  increases  linearly  with  altitude  yielding  ray  paths  which  are  arcs  of  a circle. 
The  radius  is  determined  by  the  sound  speed  gradient.  It  is  natural  then  to  break  the  "age" 
integral  into  two  parts  corresponding  to  these  two  different  ray  paths.  For  altitudes  below  the 
break  point,  only  the  first  part  of  the  integral  applies,  with  the  upper  limit  set  to  the  altitude  in 
question. 


The  density  formulas  are  given  in  the  section  describing  the  atmosphere.  The  remaining 
unknown,  the  ray  tube  area,  is  discussed  below.  The  dimensions  of  the  ray  tube  area  are  com- 
pletely arbitrary  as  they  appear  only  as  a ratio  in  all  of  the  formulas.  These  areas  may  be  found 
by  quadrature  once  the  rays  have  been  found.9  The  formula  for  rays  along  the  flight  tracks  for  a 
plane  moving  in  steady  level  flight  is  (see  Hayes  et  al9) 


sin  9„  r cos  9 
cos  0„  •'z  sin0 


d(-Z) 


(3-15) 


The  rays  are  given  in  two  segments.  The  first  has  cos  9 and  sin  9 constant.  In  the  second 
region  cos  9 = cos  0„(1  + S (Z  - Zk)]  which  is  the  Snell’s  law  solution  for  a linear  sound 
speed  profile.  Thus: 


A 


Z 


Zk  + 


sin  9„ 

6 cos2  9„ 


[sin  0„  — sin  0] 


for  Z < Zk 
for  Zk  < Z < Z, 


(3-16) 
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At  this  time,  a transformation  of  variables  is  convenient.  We  introduce: 


£ - — — A-  and  — = cos  9 = cos  9„  [1  + 8 (Z  - Zk)] 

2H  t T) 


After  considerable  algebraic  manipulation,  we  get 

T e^°  (2//,) l/ 


(3-17) 


f = 


; 


cos  « , 


c'2  p,1/2  ( sin  9„)3/1  (cos  0„)52 
sin  9„  (cos  9n)32  e 

+ JThTs 

cos  0n 

^ 3/2 


2H,8 


C°5«  (tj2  - 1) 


3/4 


2H26 

t)~  ~ e dr\ 


8 cos*  9„ 


Z\  + sin  90 


r,  - n/t,2  - 1 


TT 


(3-18) 


sin  9 „ 

The  integral  to  be  evaluated  in  the  "isothermal"  region  below  Z*  is  of  the  form 


(3-19) 


This  is  Dawson’s  integral19  which  cannot  be  further  reduced.  The  function  is  tabulated  or  can 
easily  be  evaluated  from  its  power  series 

r 


,i{°  2 \,?o  (2«  + 1)  «! 


The  integral  to  be  evaluated  above  ZA  is  of  the  form 

hi*}  = Sr, 


(3-20) 


* l /cos 
f V, 


( V 

e-Bv 

(v  + D3/4 

(aizjH 

°V 

1/2 

(v  - 1)3/4 

d-t) 


(3-21) 


where  r\  will  vary  from  1 (when  Z = Z,)  to 


1 


cos  9„ 


(when  Z = Z^).  Of  most  interest  is  the 


case  when  r/  = 1,  in  which  case  the  integrand  has  an  integrable  singularity  at  the  end  point. 
The  analysis  must  properly  account  for  this  singularity. 
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The  most  significant  contributions  to  the  integral  come  from  smaller  values  of  tj  where 
the  exponential  and  (tj  - 1)~3/4  terms  dominates.  We  rewrite  /2  as 


(t,  + 1)3/4  1 ~ 


(rj2  — l)1 


(t)2-  l)1 


G(t))  is  slowly  varying  except  near  tj  = 1 where  [1 *— — 1 2 has  an  infinite  derivative. 


To  isolate  the  effect  of  the  infinite  slope,  we  subtract  from  G(  tj)  the  quantity 
(tj  — 1 ) 1 /2/ ( 2 D 23/4).  Its  derivative  has  the  same  limit  as  that  of  <J(rj)  as  tj  — * 1.  Hence, 


(t>)  = X 


».  f(n)  e8”  . . fO)  e-8” 

dl>  + -57T  


■».  (tj  — 1 ) 324  ' 2Z)  (tj  - 1) 1/4 


77  dr,  (3-24) 


F(rj)  = <7(tj)  — 


(tj-1)1'2 


The  function  F (tj)  is  slowly  varying  over  the  entire  range  of  possible  values  for  tj,  1 < tj  ^ 2. 
We  can  therefore,  to  a good  approximation,  replace  F(r,)  by  the  straight  line 


F(r,)  = F{\)+2  F(y)-f(l)  (i>  - 1) 


(3-26) 


which  is  exact  at  tj  = 1 and  tj  = 3/2.  Equation  (3-24)  becomes 


/*  I /cos  U p-  Brt 

, zf: 


The  integrals  in  Eq.  (3-27)  can  be  evaluated  exactly  in  terms  of  the  incomplete  gamma  func- 
tions 


/ 2 <T) ,)  = 


B 1/4 


fd)  r d/4,  bv ) + 


F(l) 


+ 2lfl3/2)  - fun  r (s/4. 
B 


lDBm 

rj  = 1/cos  0O 


r (3/4,  firi) 


(3-28) 


The  incomplete  gamma  functions  T (1/4,  fir/)  and  T (3/4,  fir/)  can  be  evaluated  using  either 
the  power  series  or  the  continued  fraction  representations.  T (5/4,  fir/)  can  be  expressed  as  a 
sum  of  F (1/4,  fir/)  and  an  exponential.20  Equation  (3-28)  agrees  with  the  numerical  integra- 
tion of  Eq.  (3-21)  to  within  2%  for  all  values  of  17, . 


If  we  wish  to  consider  winds,  then  we  no  longer  get  a simple  relation  between  sound 
speed  and  cos  9.  Snell’s  law  gives  the  result  that 

cos  8 = -r1-  (3-29) 

c + u 

If  the  wind  speed  u is  a constant,  then  the  simple  law  applies.  We  can  get  a lowest  order 
approximation  to  the  effect  of  wind  by  assuming  the  wind  to  be  zero  in  the  lower  region  and 
equal  to  a constant  in  the  linear  sound  speed  region.  This  gives  a discontinuity  in  the  ray  slope 
which  violates  the  assumption  of  geometric  acoustics.  What  actually  happens  is  that  the  wave 
will  be  partially  reflected  from  this  discontinuity.  If  the  discontinuity  is  small,  then  the 
reflected  part  will  be  small.  The  solution  we  obtain,  in  any  case,  overestimates  the  signal 
strength  and  thus  gives  an  upper  bound. 


If  we  wish  to  consider  rays  other  than  along  the  flight  tracks,  the  area  formulas  must  be 
modified.  The  complication  is  not  too  severe  if  we  require  the  winds  to  be  zero  for  these  cases 
Then  the  area  formula  is  (see  Ref.  9) 
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Mach  angle.  From  the  relation 

sin  6n  — cos  n cos  <b  (3-32) 


the  angle  8„  can  be  determined  for  a given  Mach  number  and  azimuthal  angle  <b.  The  heading 
of  the  ray  with  respect  to  the  flight  tra;k  can  be  computed  from  the  relation 

tan  v = sin  d>  cot  n (3-33) 

The  general  formulas  are  given  in  Reference  9.  We  have  considered  only  the  limited 
cases  of  straight,  level,  and  steady  fight  to  simplify  the  integrals.  The  more  general  formulas 
would  require  a numerical  calculation.  Estimates  of  the  primary  effects  we  seek  can  be  found, 
however,  from  these  simpler  relations. 

Examining  the  integral  for  the  age,  we  note  that  the  biggest  variation  in  the  integrand  is 
due  to  the  exponentially  decreasing  density  with  altitude.  Thus,  the  largest  part  of  the  contri- 
bution to  the  integral  will  come  from  the  region  of  lowest  density  near  the  top  of  the  ray  path. 
Examining  Figure  3-2  which  shows  a plot  of  log  density  along  the  ray  path,  we  see  that  the 
greatest  part  of  the  path  is  in  the  region  of  lowest  density.  For  this  reason,  the  aging  of  the 
wave  takes  place  almost  entirely  near  the  top  ot  the  ray  path. 


* 
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4.  Atmospheric  model 


i 


The  ray  paths  are  determined  by  the  sound  speed  and  wind  profiles  in  the  atmosphere. 
The  altitude  at  which  the  ray  turns  over  is  determined  by  where  the  phase  velocity  (sound  plus 
wind  speed)  is  equal  to  the  aircraft  speed  (if  we  assume  no  winds  at  the  aircraft  altitude). 
Below  20  km,  data  for  the  temperature  and  wind  velocity  are  routinely  taken  by  the  weather 
service.  Sounding  rocket  data  has  provided  good,  seasonally  averaged  data  for  altitudes  up  to 
50  km.  Above  50  km  the  wind  data  is  scarce  and  many  of  the  wind  profiles  must  be  inferred 
from  models  of  the  thermosphere  and  upper  stratosphere.  In  the  regions  above  20  km,  where 
daily  weather  data  is  unavailable,  there  may  be  considerable  daily  variation  in  the  winds  due  to 
tidal  effects.  In  regions  above  100  km  the  variations  exceed  the  mean  values  so  that  no  reliable 
wind  data  for  a given  time  is  available.  Wind  speeds  of  the  order  of  100  m/s  are  found  at  125 
km  and  decrease  to  around  50  m/s  at  150  km.7  Since  no  predictive  capability  seems  reasonable 
based  on  these  wind  variations,  it  was  felt  best  to  ignore  the  wind  effect  except  to  determine 
worst  case  results  based  on  known  maximum  winds.  Otherwise,  the  results  were  computed 
ignoring  the  effect  of  winds.  Measurements  of  the  temperatures  and  density  profiles  have  been 
more  extensive  in  the  upper  atmosphere,  and  a number  of  models  exist  that  predict  tempera- 
ture and  specie  density  as  a function  of  altitude.  In  the  region  above  120  km  the  temperature 
and  density  profiles  are  strongly  influenced  by  the  solar  flux  of  energy  into  the  thermosphere. 
Variations  of  several  hundred  degrees  exist  between  the  maximum  daytime  and  nighttime  exos- 
pheric temperature.  The  minimum  exospheric  temperature,  which  is  also  used  in  determining 
the  maximum  temperature,  is  also  a function  of  solar  activity.  These  temperatures  are  also  a 
function  of  latitude  and  time  of  day.  It  is  possible,  therefore,  to  generate  a very  large  parame- 
ter range  of  temperature  and  wind  profiles  which  will  alter  the  ray  paths  and,  most  importantly. 
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the  maximum  altitude  of  the  ray.  We  have  chosen  to  find  solutions  for  a reasonable  set  of 
parameters  and  then  investigate  how  these  solutions  might  independently  vary  with  the  parame- 
ters. 


Models  exist  which  can  predict  the  temperature  and  specie  concentrations  for  a given  solar 
activity,  latitude,  and  time  of  day.  The  JACCHIA  71 21  model  is  the  one  we  are  currently 
using.  To  narrow  the  problem  somewhat,  we  have  chosen  a solar  flux  for  December  2,  1977  as 
a representative  value.  This  also  corresponds  to  one  of  the  "event”  days.  We  chose  a latitude 
of  41°N  and  a time  of  12  noon  to  correspond  to  Concorde  arrival  in  New  York. 

The  important  parameter  that  the  models  predict  in  terms  of  the  ray  paths  is  the  local 
sound  speed  rather  than  the  temperature.  At  altitudes  above  100  km  the  composition  of  the 
atmosphere  changes  resulting  in  a lower  mean  molecular  weight.  This  has  the  effect  of  raising 

the  sound  speed  for  a given  temperature  since  c2  = This  can  <ncrease  *he  sound  speed 

M in 

by  as  much  as  20%  more  than  that  obtained  using  a constant  mean  molecular  weight  atmo- 
sphere at  180  km.  This  effect  lowers  the  peak  altitude  of  the  ray  by  20  km. 

In  Figure  4-1  we  show  a comparison  of  the  JACCHIA  ’71  mode!  sound  speed  compared 
to  the  U.S.  1976  Standard  Atmosphere  profile.22  The  results  are  similar  enough  that  either 
could  be  used.  We  use  the  JACCHIA  ’71  model  results  for  the  calculations  presented  here. 

In  order  to  make  our  calculations  analytically  integrable,  we  have  chosen  to  make  another 
simplification  of  the  atmospheric  model.  We  construct  a model  atmosphere  that  has  a constant 
sound  speed  up  to  some  turning  point  Z;  and  then  a linear  increase  in  the  sound  speed  up  to 
the  altitude  Z,  where  the  ray  reaches  its  maximum  altitude.  The  constant  sound  speed  is  equal 
to  the  sound  speed  c„(295  m/s)  at  the  Concorde  cruise  altitude.  The  point  ZA  (103  km)  is  the 
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Fig  4 1 - Comparison  of  ihe  sound  speed  from  Ihe  Jacchi.i  model  atmosphere 
for  a midday,  winter,  midlatitude  period  of  low  solar  activity  with  the  1976  U S 
Standard  Atmosphere. 


highest  altitude  at  which  c„  occurs  on  the  actual  sound  speed  profile  as  shown  in  Fig.  4-2. 
From  Zk  to  Zr,  the  turnaround  altitude,  the  sound  speed  is  assumed  to  increase  linearly  from 
295  m/sec  to  590  m/sec.  That  is,  sound  speed  is  modeled  by 

c = c„  0 ^ Z < Zk  (4-1) 

c - c0  (I  + S(Z  - Zk)  1 ZA  < Z < Z, 


(1  - cos  9o) 
cos  9„(Zi  - Zk) 


(4-2) 


From  Fig.  4-2  we  see  that  the  assumed  profile  is  in  good  qualitative  agreement  with  the  actual 
profile.  By  using  this  simplified  model,  we  reduce  the  ray  paths  to  straight  lines  below  Zk  and 
arcs  of  circles  above  Zk.  The  actual  sound  speed  profile  above  Zk  is  not  well  known,  and  since 
our  assumed  straight  line  fits  the  best  guess  rather  well,  we  will  not  concern  ourselves  with  the 
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SOUND  SPEED  (m/sec) 

Fig  4 2 — Comparison  of  sound  speed  for  besl  guess  aimosphcre 
wilh  ihe  assumed  sound  speed  profile 

accuracy  of  our  model  atmosphere  above  Z — Zk.  To  lest  ihe  validity  of  our  model  below 
Z “ ZA  we  performed  the  ray  tracing  for  the  U S.  Standard  Atmosphere,  for  several  initial 
angles.  The  results  are  shown  in  Figs.  4-3  and  4-4.  For  = 60°,  which  corresponds  to  the 
on-axis  case  for  a Mach  number  of  2,  we  see  that  our  assumed  straight  line  ray  differs 
insignificantly  from  the  actual  ray.  This  is  also  the  case  for  ft,,  = 50°.  For  ft,,  = 40“  the  length 
of  the  two  rays  is  nearly  the  same  so  the  computed  age  and  hence  the  final  amplitude  and  fre- 
quency will  be  nearly  the  same.  This  is  especially  true  since  the  vast  majority  of  the  aging 
occurs  above  Z*.  We  will,  however,  in  this  case  underestimate  the  propagation  range  by  10 
km.  For  <t„  = 30°  we  would  probably  have  some  difference  in  the  age  and  a significant  error 
(60  km)  in  the  propagation  distance.  An  enormous  error  occurs  when  ft.,  is  less  than  26°  since 
in  this  case  the  wave  never  reaches  the  thermosphere  but  is  reflected  from  below  50  km.  Such 
a ray,  however,  would  not  reach  the  ground.  Thus,  with  (he  U.S.  Standard  Atmosphere,  little 
error,  except  for  range  predictions  for  ft  ^40°,  is  incurred  by  using  an  isovelocity  profile  below 
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Z When  Winds  blowing  in  the  direction  of  the  aircraft  flight  are  present  in  the  ozonosphere, 
for  cases  far  enough  off  track,  we  may  predict  thermospheric  returns  when  the  return  may 
actually  be  coming  from  below  60  km. 

The  pressure  and  density  for  our  model  are  found  by  integrating  the  equations  for  hydro- 
static balance. 

<«> 

dZ  c- 

Again  we  break  the  integral  into  two  regions,  one  where  the  sound  speed  is  constant,  and  the 
other  where  c = i,.ll  + 6 (Z  — Z,)J.  Then, 


-1— 1 
pip  = e 1 1 

for  Z < Z< 

Z-ZK 

M, ll  -*1/  - A'l 
PIPk  = f ' 

for  Z;  < Z < Z 

where  H,  and  H:  are  scale  heights.  The  density  is  found  from  the  relation  P - Thus 


P = P e 


Z < ZK 


(4-4) 


2 - \ 

" «,[1  -*'Z-Z„>) 

P = P'  [1  + 6<Z  - Z.)]y  Z' 
Since  the  scale  height  is  not  constant 


< Z < Z pt  - P.  c 

with  altitude,  using  the  scale  height  computed  from 


H i „ 21  will  lead  to  an  error  in  p and  p at  the  final  altitude.  To  overcome  this  we  compute 

c,r 


two  scale  heights  H]  and  H:  which  are  computed  by  matching  the  pressures  at  the  turning  point 
Z.  and  the  final  altitude  Z..  This  has  the  same  effect  as  redefining  an  effective  g for  each 

region. 


In  the  case  with  constant  winds  we  have  from  Snell  s law  that 


(c  u 1 cos  f > - c - c (1+  8(Z  - ZK)  1 


(4-5) 
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Solving  for  the  top  of  the  ray,  where  cos  9 *=  1 , we  can  find  the  sound  speed  lapse  rate  8.  Both 
the  model  atmospheric  properties  and  the  ray  equation  are  now  completely  defined.  From 
Snell's  law 


cos  9 - 


cn[\  + 8 (Z  - Zk)  1 
(c'  + u) 


(4-6) 


With  the  wind  equal  to  zero  at  the  Concorde  fight  altitude,  we  have 


cos  0„ 


11 


r- 


I 


5.  Linear  Attenuation 


Acoustic  attenuation  in  air  at  infrasonic  frequencies  is  negligible  at  sea  level  (<  .01 
db/km).23  At  high  altitudes,  however,  the  attenuation  of  infrasound  can  become  quite  large 
because  the  acoustic  attenuation  per  wave  length  is  a function  of  ftp24  Thus  at  160  km,  where 
pip K = 3 x 10-*.  the  attenuation  per  wavelength  at  1 Hz  is  equivalent  to  the  sea  level  attenua- 
tion per  wavelength  at  333  MHz.  (The  subscript  g will  refer  to  values  at  sea  level.)  There  are 
four  major  sources  of  acoustic  attenuation  in  air;  viscosity,  heat  conduction  and  the  rotational 
and  vibrational  relaxation  of  02  and  N2-  In  this  section  we  consider  each  of  these  mechanisms 
and  obtain  an  attenuation  model  for  the  upper  atmosphere. 

A.  Vibrational  Relaxation 

Vibrational  relaxation,  which  is  the  dominant  attenuation  mechanism  at  low  frequencies  at 
atmospheric  pressure,  is  extremely  sensitive  to  the  presence  of  trace  constituents  in  the  air, 
mainly  water  vapor  and  C02.  For  example,  the  attenuation  at  1 kHz  at  10%  relative  humidity  is 
twice  that  at  20%  relative  humidity  while  at  20  kHz  this  relation  is  reversed.23  We  shall  show, 
however,  that  for  the  effective  frequencies  (f/p)  of  interest,  these  complex  vibrational  relaxa- 
tion processes  are  frozen  out. 

Of  primary  concern  is  the  vibrational  relaxation  of  02  since  the  relaxation  frequency  of  N2 
is  orders  of  magnitude  lower.25  At  high  humidities  the  vibrational  relaxation  frequency  may 
exceed  100  kHz.  However,  above  100  km  the  concentration  of  water  vapor  is  much  less  than  1 
ppm22  so  we  are  interested  in  the  relaxation  frequency  for  dry  air  which  is  known  to  be  between 
30  and  40  Hz/atm.23  2t  In  the  absence  of  water  vapor  the  vibrational  relaxation  frequency  is 
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strongly  influenced  by  the  presence  of  minute  quantities  of  C02.  For  C02-free  air  the  relaxa- 
tion frequency  is  6 Hz/atm.  Though  the  mixing  ratio  of  C02  at  20  km  is  320  ppm22,  it  is  rea- 
sonable to  expect  it  to  be  significantly  lower  above  100  km.  Thus,  we  expect  the  relaxation  fre- 
quency above  100  km  to  be  6 Hz/atm,  though  it  could  conceivably  be  as  high  as  30  Hz/atm 
(which  corresponds  to  300  ppm  C02).  The  lowest  frequencies  that  we  will  be  concerned  with 
are  of  the  order  of  .01  Hz  which  corresponds  to  30  kHz/atm  at  100  km  and  3 MHz/atm  at  150 
km.  Since  these  frequencies  are  many  orders  of  magnitude  larger  than  the  largest  conceivable 
value  for  the  relaxation  frequency,  we  can  safely  conclude  that  vibrational  relaxation  can  be 
neglected. 


B.  Classical  Attenuation 


Heat  conduction  and  viscosity  are  the  dominant  attenuation  mechanisms  in  the  thermo- 
sphere at  infrasonic  frequency.  Together  they  comprise  what  is  referred  to  as  classical  attenua- 
tion. The  classical  attenuation  coefficient  for  air  is  given  by 


r2 


8 


f2  3pcJ 


1,4 

y - 1 

K 

4 

c„ 

(5-1) 


where  tj  is  the  shear  viscosity,  Cp  the  constant  pressure  heat  capacity  and  k the  thermal  conduc- 
tivity. Since  p = pc2/y , we  can  rewrite  Eq.  (5-1)  as 


“,/P  = 

f2  3yc 


V + 


l 


C„ 


(5-2) 


Thus  ad  is  directly  proportional  to  the  square  of  the  frequency  and  inversely  proportional  to 
ambient  pressure.  Equation  (5-2)  is  valid  below  the  translational  relaxation  frequency,  i.e.  as 
long  as  the  wavelength  is  greater  than  a mean  free  path.  From  the  kinetic  theory  of  gases,  the 
thermal  conductivity  and  viscosity  are  related  by  the  Euchen  relation.27 


9y  - 5 


4y 


r 19  _ 
CnV  - C„V 


(5-3) 
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Thus  from  Eqs.  (5-2)  and  (5-3) 


a./P  394  7T3 


f2  l 35  > J c 

so  that  the  classical  attenuation  depends  solely  on  the  ratio  tj/c.  The  temperature  dependence 
of  the  viscosity  is  given  by  Southerland’s  formula 

_ b<Jm 
V T -I-  110.4 

where  b„  = 1.458  x 10~6  kg/(sec-meter-°K12).32  If  is  the  viscosity  at  sea  level  then  at  160 


V = 1.856  ifK 


so  that 


c = 1.735  ce 


* = 1.07  -2! 
c c„ 


from  which  we  conclude  that  a^p/f2  is,  to  a good  approximation,  independent  of  temperature 
and  hence  altitude.  We  shall  therefore  use  measured  values  of  ac\p/f 2 appropriate  to  20°C. 


C.  Rotational  Relaxation 


The  rotational  relaxations  in  02,  \'2  and  air  behave  like  a simple  single  relaxation  process. 
For  such  a relaxation  process  the  attenuation  coefficient  is  given  by 


2(oi  mlK ) r 


a.i  = 


where  is  the  rotational  relaxation  frequency  at  one  atmosphere  and  (ar„,X)m„  is  the  max- 


imum value  of  a„„K,  which  is  attained  when  f = frm.  When  the  effective  frequency,  — -,  is 

P 

less  than  frnr  Eq.  (5-6)  becomes 


30 


(5-7) 


**  rot  P_  2(  & rot  ^ m dxP/f 

f 2 f„„C 

Comparing  Eqs.  (5-2)  and  (5-7),  we  see  that  a,,  and  am,  both  have  the  same  dependence  on 
pressure  and  frequency. 

The  temperature  dependence  of  the  right  hand  side  of  Eq.  (5-7)  presents  some  difficulties. 
It  is  known’1*  that  the  rotational  relaxation  frequency,  for  N,  is  temperature  dependent 
decreasing  from  over  300  MU/  at  295IIK  to  around  75  MHz  at  773°K.2i<  This,  however,  is  par- 
tially offset  by  the  increase  of  sound  speed  and  the  decrease  of  («„((x)mjx  with  increasing  tern 
perature.  Since  the  contribution  of  to  the  total  attenuation  is  relatively  small,  and  the  data 
at  high  temperatures  is  sparse,  we  elected  to  ignore  the  temperature  dependence  of  It  can 
be  shown  that  this  will  result  in  underestimating  it  by  no  more  than  25%  at  160  km  and  by  con- 
siderably less  ul  lower  altitudes. 

D.  Empirical  Value  for  a 

With  the  vibrational  relaxation  of  02  frozen  out,  the  attenuation  in  air  and  N2  are  virtually 
identical.  Since  the  air  below  160  km  is  primarily  N2,  we  shall  use  a value  for  a appropriate  to 
a nitrogen  atmosphere  In  Fig.  5-1,  taken  from  Greenspan,27  experimental  and  theoretical 
values  for  a\  for  N2  are  plotted  as  a function  of  Reynolds  number.  The  values  for  the  Rey- 
nolds number  appropriate  to  our  problem  at  i 60  km  and  120  km  are  indicated  by  the  arrows.  It 
is  evident  that  we  are  operating  well  within  the  region  where  ap/f 2 is  a constant.  That  is,  we 
are  well  above  the  Reynolds  number  at  which  either  rotational  or  translational  relaxation 
occurs.  (With  regard  to  the  latter,  at  160  km  the  mean  free  path  is  53  m while  our  pulse  length 
is  about  4 km.) 

The  values  quoted  in  the  literature  for  a//2  for  N,  at  atmospheric  pressure  are  all  in  rea- 
sonably good  agreement  with  one  another.  Values  of  1.79  x.  10  11  Neper  • sec:/m  (Sessler),2g 
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Fig.  5.1  - Attenuation  per  unii  wavelength  and  sound  speed  for  molecular  ni- 
trogen are  shown  as  a lunction  of  Reynolds  number  Reynolds  numbers  ap- 
propriate to  the  thermosphere  for  our  problem  are  shown  by  the  arrows  From 
Greenspan,  J Acoust.  Soc.  Am  31.  Id5(l 


1.66  x 10  " Neper  • secVm  (Winter  and  Hill),2"  1.64  x 10"  Neper  ■ sec-’/m  (Parker  el.  at.)* 
and  1.74  x 10  " Neper  sec  7 m (Greenspan)-’7  were  averaged  to  obtain 


a = 1.70  x 10 


r Neper/m 


(5-8) 


wtth  a standard  deviation  of  .07.  This  is  the  value  for  the  attenuation  coefficient  that  we  will 
use  at  all  altitudes. 
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6.  Linear  Solution  and  Crossover  Model 


When  the  acoustic  attenuation  coefficient  is  proportional  to  /2,  as  in  Eq.  (5-8),  the  Blo- 
khintsev  variable,  in  the  absence  of  winds,  satisfies  the  following  reduced  wave  equation  along 
a ray 


dK+  ry  ty_a{s)  djv 

ds  c2JpcA„  Bt  dr1 


(6-1) 


Here  s is  the  distance  along  the  ray,  r is  the  retarded  tirnc,  and  the  attenuation  factor  o(s)  is 
defined  by 

a(s) 


a(s)  = 


(6-2) 


The  attenuation  factor  depends  on  s due  to  the  pressure  dependence  of  a.  Since  the  particle 
velocity  is  related  to  the  Blokhintsev  variable  by 

y 


q = 


y/pcA,, 


Eq.  (6-1)  can  be  rewritten  as 


(6-3) 


(6-4) 


For  an  N wave,  where  the  time  of  the  positive  phase  is  given  by  the  relative  nr.por- 


r q dv 


b2V 


tance  of  the  nonlinear  term  and  the  linear  attenuation  term  a(s)  — j is  measured  by 

cl  OT  OT 


the  ratio 


R = 


Yqtp 

ac 2 


(6-5) 


which  is  called  the  Gol’dberg  number.  When  R is  very  large,  the  linear  attenuation  term  can 
be  neglected,  yielding  the  inviscid  Burgers'  equation 
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9s  c2Vp^7  8t 

which  was  solved  in  Section  3.  When  R is  small,  the  nonlinear  term  can  be  neglected  and  the 
propagation  is  governed  by  the  diffusion  equation 

9 ^ , d2Y  ,, 

= a (s)  —7  (6-7) 


We  assume  that,  at  a given  altitude,  either  Burgers’  equation,  Eq.  (6-6),  or  the  diffusion  equa- 
tion, Eq.  (6-7),  applies  This  is  possible  because,  due  to  the  increasing  a and  c,  the  Gol’dberg 
number  decreases  rapidly  with  altitude  so  the  region  where  the  complete  Eq.  (6-4)  is  required  is 
very  narrow.  An  appropriate  place  to  make  the  transition  from  Eq.  (6-6)  to  Eq.  (6-7)  is  where 
the  incremental  decrease  in  energy  due  to  nonlinear  attenuation  becomes  equal  to  the  incre- 
mental decrease  in  energy  due  to  linear  attenuation.  That  is,  we  switch  from  Eq.  (6-6)  to  Eq. 
(6-7)  when 


± M 

E 9s 


1 M 

E 9s 

' NONLINEAR 


The  energy  in  the  N wave  is  proportional  to  qptn>  where  qr  is  the  peak  value  of  the  velocity  and 
ir  is  the  time  of  the  positive  phase.  From  the  N wave  solution  we  obtain  at  any  point  s 


1 M 

E 9s 


-rg , 

NONLINEAR  2 C2ln 


The  spectrum  of  the  A wave  has  its  maximum  at  fp  ~ — — , so  the  linear  energy  attenuation 


coefficient  is  approximated  by 


— 2 a = — 8 it  2af} 


(6-10) 


Equation  (6-8)  is  thus  satisfied  when 


(6-11) 


We  shall  therefore  use  Eq.  (6-6)  when  R > 17.5  and  Eq.  (6-7)  when  R < 17.5.  [We  note 
that  Eq.  (6-4)  can  be  solved  exactly  when  a is  independent  of  s.  We  have  examined  the  solu- 
tion and  found  it  consistent  with  this  choice  for  the  cross  over  value  for  R.  At  R - 30  the 
wave  is  definitely  a weak  shock  and  at  /?  =>=10  it  is  definitely  linear.]  We  shall  examine  the 
consequences  of  modeling  Eq.  (6-4)  in  this  way  later  on  in  greater  detail. 


Since  the  weak  shock  solution  for  an  N wave  is  always  an  N wave,  the  initial  condition  for 
our  solution  to  Eq.  (6-7)  will  be  an  N wave.  We  assume  that  at  the  transition  point  sL,  the 
wave  form  is  given  by 


V = - 


- rL  < T <Tt. 


(6-12) 


If  we  change  variables  from  s to 


j"  a(s)ds 


(6-13) 


then  Eq.  (6-7)  becomes 


0 V __  tfv 

9*  dr2 

which  is  a simple  diffusion  equation  whose  solution,  subject  to  the  given  initial  condition,  is 


(6-14) 


1 * 

erf 

r + 1 
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T - 1 

* 2 T| 

— erj 

2fe 

_ ±L 

T + 1 

■Jir 

exp 

2 U 

- exp 


- 1 
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where  V - V/VL,  f - r/rL  and  the  "linear  age"  is  defined  by 


(6-15) 


(6-16) 


The  waveform  is  plotted  for  several  values  of  in  Fig.  6-1.  For  > 1.  V approaches  the 
form 


f -rJ/4£^ 
6-v/ir^f  P 


(6-17) 
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Fig  6 1 - Waveform  vs  retarded  time  for  several  values  of 
the  age  parameter  showing  effect  of  linear  aging. 


' 


Equation  (6-17)  has  a maximum  at  f,„  = 7Zf  so 

Similarly,  the  energy  is  given  by 


and  the  peak  frequency  by 


E _ /T  1 
El  “ V rr  24f| 


U > 1 


(6-18) 


(6-19) 
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peak 


l 


(E  ^ I 


(6-20) 


ft  n/2££ 

From  these  asymptotic  forms  and  numerical  analysis  of  Eq.  (6-15),  we  can  obtain  simple  empir- 


ical formulas  for  and  E/EL  as  a function  of  f £ 


peak 


El 


(1 


1 f ,452c  l4*f 
[1.452  - y-  +6 

1 - .17^/rf  Sff 

+ 2.5ff  + 1.3«,2  + 24 


(6-21) 


(6-22) 
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(6-23) 


where  t - 2.0815760  and  /,  - — 1 — . Equations  (6-21)-(6-23),  which  are  plotted  in  Fig.  6-2, 

are  accurate  to  within  1%  for  all  values  of  f£.  All  of  the  pertinent  information  about  the  linear 
solution  can  easily  be  obtained  from  Eqs.  (6-14),  (6-20),  (6-21),  and  (6-23)  once  one  deter- 
mines the  linear  age  f£.  The  problem  is  thus  solved  once  one  obtains  \ fr°m  Eq  (6-13). 


C, 


Fig  6 2 — Decay  of  amplitude,  peak  frequency, 
and  wave  energy  with  increasing  linear  age 


The  dependence  of  the  attenuation  factor  on  s is  due  solely  to  the  dependence  of  the 


ambient  pressure  p on  s.  Hence 


a(s) 


aKPn 

p(s) 


(6-24) 


where  aK  and  pK  are  the  attenuation  coefficient  and  pressure  at  ground  level. 


In  order  to  determine  the  value  of  the  Blokhintsev  variable  V at  the  turning  point,  we 
evaluate  *„  which  is  obtained  from  Eq.  (6-13),  by  integrating  a(s)  along  the  ray  from  the  alti- 
tude Z,  (where  R = 17.5)  to  Z,  (the  turning  point  altitude).  It  turns  out  that  ZL  will  always 
be  above  the  "knee"  in  the  sound  speed  profile,  which  is  located  at  Zk  = 103  km.  From  our 
previous  analysis,  the  static  pressure  is  given  by 


p(Z ) = p(Zk)  exp 


-(Z  - Z*)cos  90 

h2 

cos  9„  + (1  - cos  9„) 

N 

1 

N 

Z,  - Z*  j 

(6-25) 


where  the  scale  height  H2  is  defined  by 


H, 


(Zf  - Zk) 


2 In  [p(Zk)lp(Zf)\ 

The  ray  is  an  arc  of  a circle  of  radius  ( Z f - Z*)/(  1 - cos  0„)  so  that 


(6-26) 
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(6-27) 


p(d)  = p(Zk)  e 


IZ,  - Zt)  cosW„  i _ cos*),, 
H2(\  - cos0o)  cost* 


The  limits  of  integration  are  6 — 0 and 


(6-28) 


9 i = cos  1 


cos0„  + (1  - COS0„) 


ZL  ~ Zk) 

z,  - zk 


From  Eqs.  (6-13),  (6-24),  (6-27)  and  (6-28)  we  obtain 


Z,  ~ zk 

1 - cos#,, 


a*p« 

p(Zk) 


Bqos90 

e cos"  d6 
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(6-29) 


(6-30) 


where 


B 


( Z,  - Zk)  cos#,, 


We  now  change  variables  from  # to  £ 


H2(  1 - cos#,,) 

1 


(6-31) 


cos# 


- 1 which  gives 


X, 


axPn  (2/  — Zk)  8(  l-cos#,,) 
e 


J.  COS»L 
0 


p(Zk)  [1  - cos#,,] 


dl 


(4  + 1)  (4  + 2) 1/2  4 


1/2  y 1/2 


(6-32) 


From  Eq.  (6-28) 


eB(i-cos»„)  = p(zk)/p{Zf) 


(6-33) 


and  since  aKpK/p(.Z,)  = a(Z,),  Equation  (6-32)  becomes 


I -cos#, 


Xi 


a(Z,) 


/•  COS  H. 

(z,-z,i  r 


-8cos#„{ 


rf4 


(6-34) 


[l-cos#„]v  / Jo  (4  + 1)  (4  + 2),/2  4,/2 

Since  B cos  #„  is  a fairly  large  number,  the  most  important  contributions  to  the  integral  come 


from  small  values  of  4.  We  can  thus  make  the  substitution 


(4  + 1)  (2  + 4) 1/2  = n/2 


1 + 4"  4 


>/2  e5</4 


(6-35) 


and  obtain 


1 -8tOS«„+f  { 

_ a(Z,)  (Zf~  zk)  r cos*,  - ,1 

*'  “ [ 1 — COS#,,]  \/2  Jo 


4 


1/2 


(6-36) 


or 


Xt 


a(Zf)  ( Zf-Zk ) 


[1  - cos#,,] 


flcos#,,  + — 


The  Blokhintsev  variable  at  Z = Z,  is  given  by  Eq. 


16-37) 


Having  obtained  Eq.  (6-37),  we  are  now  ready  to  reexamine  our  treatment  of  the  transi- 
tion from  weak  shock  theory  to  linear  theory.  Due  to  the  properties  of  the  error  function, 
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lu,„s  on.  be  nearly  independen,  .r  Thar  is.  unless  «r  Is  very  nearly  equal  .0  aero,  .he 
error  function  is  almos.  exactly  equal  10  unity.  Hence,  in  general. 


a(Z,)  (Z,  - Zo 
Xr  ~ (1  - cos0„) 


2\B  cos  0„  + 


(6-38) 


The  same. .inn  elfec,  comes  .hour  because,  due  .0  .be  na.ure  of  .he  -ay  - — — 
mos,  of  us  „me  a.  the  hrghe,  al-Uudes  .here  .be  a.r.nua.ron  is  mos,  severe.  Thus,  mos,  of  .he 
aging'  occurs  near  .be  ,00  of  .be  oa.h  and  hence.  ,.  is  rela.ively  unimoor.an,  e.ac.iy 
„„e,e  ,he  ...e  s.aned.  As  a resul,  of  .bis  sarura.ion  effect.  x,  - 

independent  of  .ba,  .«  choo*  for  .he  crossover  va.ue  for  ,he  Co.dberg  number.  From  K 
<6-,g)  the  peak  value  of  .be  Blokbin.sev  .enable  a.  .be  .u.nm,  point  is  given  by 


rrjw_ 

i>peak(Z/)  ” V ne  6 |*f(Z,) 


or,  with  x i constant, 


T 


V^(Zf)  - const.  X [,,(Zj]  tlZi) 


(6-39) 


(6-40) 


7 is  aoverned  by  weak  shock  theory,  the  impulse,  given  by 
Since  the  propagation  up  to  Z,  is  governed  oy 


tr[Zt)  y(Z,),  is  a constant  independent  of  ZL  Hence, 

V^Z,)  “ Const,  x ,,iZL) 


(6-41) 


..  f the  nulse  t will  increase  with  increasing  ZL.  Thus,  using  a 
Due  to  the  nonlinear  stretching  of  the  pulse,  tf  win 

smaller  value  for  .be  crossover  Goldberg  number  .ou.d  resul,  m a large,  value  fo,  f 
However,  a,  ,he  wave  propaga.es,  small  changes  in  ...Unde  produce  large  changes  m ,he 
Oordberg  number  bu,  only  rela.rvely  small  changes  in  ,,  Hence.  ,be  exac,  value  .e  choose 
f„,  me  crossover  Gol'dberg  number  is  no,  crucial  Doubling  »,  halving  .he  crossover 

Oordberg  number  changes  our  result  by  less  than  2 dB 
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The  preceding  results  help  illustrate  an  important  point.  Weak  shock  waves  propagate 


much  belter  in  the  upper  atmosphere  than  acoustic  waves  of  the  same  frequency.  This  is  fairly 
evident  from  the  fact  that  the  energy  of  a weak  shock  wave  attenuates  asymptotically  like  v 1 ’ 
while  the  energy  of  a linear  <V  wave  attenuates  like  \ 1 2.  Once  linear  attenuation  takes  over, 
the  wave  attenuates  much  more  rapidly  than  it  would  have  had  the  wave  remained  shocked.  It 
is  therefore  reasonable  to  expect  that  once  the  wave  becomes  linear,  it  will  never  reshock.  We 
therefore  assume  that  once  the  wave  becomes  linear,  it  remains  linear.  Moreover,  if  we  use 
weak  shock  theory  to  solve  the  entire  problem,  we  will  obtain  an  upper  bound  for  the  solution. 

To  evaluate  Fat  some  height  Z other  than  Z,,  we  replace  the  upper  limit  in  Eq.  (6-30)  by 
the  appropriate  9 which  is  calculated  by  replacing  ZA  by  Z in  Eq.  (6-29).  For  the  downgoing 
part  of  the  ray,  we  replace  9 as  calculated  above  by  —9.  The  integral  is  evaluated  in  the  same 
manner  as 
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7.  The  Caustic 


I.  Introduction 

The  reflection  of  a wave  by  a continuously  stratified  medium  produces  a caustic  For  the 
sonic  boom  the  caustic  coincides  with  the  turning  point  for  <t>  — 0 and  occurs  somewhat  beyond 
the  turning  point  for  <t>  r*  0.  At  the  caustic  the  ray  tube  area  A„  — 0,  so  the  ray  theory 
predicts  an  infinite  pressure.  The  singularity,  however,  is  integrable  so  the  infinite  value  does 
not  affect  the  calculation  of  either  the  nonlinear  or  linear  age  parameter.  The  linear  problem 
can  be  solved  exactly,  at  and  in  the  vicinity  of  the  caustic,  in  terms  of  Airy  functions'31,  but  the 
non  linear  problem,  to  our  knowledge,  remains  unsolved.  One  would  expect  that  in  the 
immediate  vicinity  of  the  caustic,  the  linear  and  nonlinear  solutions  would  be  quite  different 
(e  g.,  the  formation  of  a Mach  stem,  self  refraction,  etc.),  so  there  would  appear  to  be  little 
advantage  to  using  the  Airy  function  solution  in  preference  to  the  ray  theory. 

We  shall,  therefore,  treat  the  caustic  by  combining  the  weak  shock  theory  with  the 
approach  of  Keller32  33  and  others34  35  3*  Keller  handles  the  caustic  by  assuming  that  when  a 
wave  passes  through  a caustic,  the  net  effect  is  to  introduce  a n/2  phase  shift  to  the  waveform 
We  therefore  compute  the  age  by  assuming  that  the  amplitude  is  governed  by  ray  theory  (since 
the  integral  converges)  and  determine  the  waveform  by  introducing  a ir/2  phase  shift  at  the 
caustic  Since  the  ray  theory  gives  an  infinite  result  at  the  caustic,  we  cannot  use  ray  theory  to 
determine  whether  or  not  the  Mach  number  in  the  vicinity  of  the  caustic  becomes  too  large  for 
weak  shock  theory  to  be  valid  To  determine  whether  the  weak  shock  theory  is  still  appropri- 
ate. we  rely  on  the  empirical  observation  that  the  focus  factor,  in  practice,  never  exceeds  five  in 
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the  case  of  a simple  caustic37  ’8.  If  the  Mach  number  is  still  sufficiently  small  when  the  focus 
factor  equals  five,  we  shall  assume  that  the  weak  shock  theory  can  be  used  through  the  caustic. 

II.  Analysis  of  "U"  Wave 

We  model  the  caustic  by  assuming  that  the  N wave  undergoes  a ir/2  phase  shift  at  the 
caustic  and  then  continues  aging.  There  is  considerable  empirical  evidence  that  this,  indeed,  is 
what  happens.  In  their  study  of  focused  booms.  Wanner  et  al37  show  many  sonic  boom 
waveforms.  In  every  case  where  a waveform  can  be  identified  with  a ray  which  has  passed 
through  a caustic,  the  waveform  has  the  U shape  characteristic  of  an  N wave  with  a ir/2  phase 
shift 


To  determine  the  effect  of  the  caustic,  we  start  with  an  N wave  with  amplitude  Vc  and 
period  2 tc  as  shown  in  Fig.  (7-1) 


V(r)  = - 

T , 


— T ^ T ^ 


The  Fourier  integral  of  Fis  given  by 


V ..  = -- 


2K 

sin((orl.) 

COs(a>T, ) 

Tc 

Jr 

CD 

(7-1) 


(7-2) 


After  passing  through  the  caustic,  the  waveform  undergoes  a W2  phase  shift  and  F(r) 
becomes  U(t)  given  by 


U (t)  — — f cos(<o  r)  V dw 
7 r Jo 


U(t)  = - — V, 

TT 


i - ~r  loel 


T + T, 


(7-3) 

(7-4) 


U(t)  and  V(t)  are  plotted  in  Fig.  7-1. 


The  function  f/(r)  has  two  real  roots  located  at  t = ±t,  with  r,  given  by 

— = .833556500  (7-5) 

Tc 
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Fig.  7.1  — U wave  results  from  a W 2 phase 
shift  of  the  shown  N wave 


When  the  "U"  wave  ages  into  a "N”  wave,  the  impulse  from  r = -»tor  = -T<  will  be  con- 
tained in  the  first  positive  part  of  the  N wave.  Let  W(t)  be  the  indefinite  integral  of  U(t). 
Then 


H/(r)  = / U( t)  (It  = 


T 


(7-6) 


The  first  positive  impulse  of  the  "U”  wave,  /y,  is  given  by 

V T2 

lv  = W (—Tr)  - H^t-oo)  = -L-L  = .3819  V'T, 

TTT, 

while  the  positive  impulse  for  the  original  N wave  is  given  by 


so  that 


V(t)  dr  = 


lv  - .76373914  /v 


(7-7) 


(7-7a) 


(1-lb) 


«*  • 
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A.  Exact  Solution 


A solution  to  the  inviscid  Burgers’  equation 

M = t/M  (7-8) 

Of  dr 

where  f '•>  the  age  and  t the  retarded  time  is 

F = F(t  + f K)  (7  9) 

For  the  U wave  at  f = 0,  F(t)  is  given  by  Eq.  (7-1).  A "IT  wave  is  shown  for  f = 0 and 

f — f 4 at  the  top  and  bottom  of  Fig.  7-2.  For  f = £4  the  wave  has  distorted  according  to  Eq. 
(7-9)  and  a shock  has  formed.  (The  waveform  will  continue  to  age  and  will  eventually  turn 
back  into  an  N wave.)  V4  is  the  peak  amplitude  and  - VB  the  shock  amplitude  when 
£ = in  the  initial  waveform,  V A occurred  at  r = —t4  and  VH  occurred  at  r = -rs.  In  the 
aged  waveform,  V , occurs  at  r = — r 4 — f 4 V 4. 

The  area  in  region  I in  the  top  of  Fig.  7-2  is  equal  to  the  area  in  region  1 in  the  bottom  of 
Fig.  7-2.  The  same  is  true  for  region  II  From  Eq.  (7-6)  t'.ien 

A , = H/,(-t,)  - W(-r  4) 

An  = W(-tb) 

While  from  Fig.  7-2,  we  can  see  that 

a,„-  V}\  l a 

Landau’s  law  or  equal  areas"1  requires  that  the  shock  be  located  so  that 

A / + An  + Am  — W(—Tr)  = lu 
so 

m-rB)  - IF ( r 4)  = y f - yl] 
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(7-10) 

(7-11) 


(7-12) 


(7-13) 


(7-14) 


£=0 


Fig.  7.2  — Nonlinear  evolulion  of  V wave 
towards  the  IV  wave  final  form 


By  definition,  when  £ = £4  the  retarded  time  for  equals  the  retarded  time  for  VB 

which  is  also  the  retarded  time  for  the  shock.  Hence, 

TB+V£B-TA+VAtA  (M5) 

Equations  (7-14)  and  (7-15)  along  with  Eq.  (7-6)  for  W and  Eq.  (7-4)  for  i can  he  solved 

numerically  to  give  the  peak  value  for  U as  a function  of  f.  The  result  is  shown  in  Table  1. 

Note  that  the  peak  amplitude  drops  to  reasonable  values  very  quickly.  For  example, 

U k = 4 vr  by  the  time  the  age  variable  £ = .000015  and  is  down  to  21/,  when  £ — 008. 

This  implies  that  we  need  not  be  concerned  about  the  logarithmic  singularity  in  the  original  L 

wave  since  it  disappears  very  quickly.  In  any  case  the  logarithmic  singularity  which  is  due  to 
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Table  I 


£ 

Exact 

Eq.  7-20 

Eq.  7-29 

001.48  x 10"5 

4.0 

1.38 

1.07 

000.0081 

2.0 

1.36 

1.06 

000.2231 

1.0 

1.11 

0.883 

001.705 

0.5 

0.60 

0.517 

009.814 

0.25 

0.27 

0.2516 

044.198 

0.125 

0.131 

0.1253 

186.35 

0.0625 

0.640 

0.06254 

the  zero  rise  time  of  the  N wave,  would  not  actually  be  present  for  a real  sonic  boom  which  has 
finite  rise  time. 


B.  Reduction  to  an  Equivalent  "N"  Wave 


The  peak  value  of  the  Blokhintsev  variable  ffor  an  N wave  is  given  by 


V = 

r peak 


K re 


(7-16) 


Te 


+ t 


For  an  N wave  whose  initial  amplitude  was  Vc  and  whose  initial  period  was  2 tc,  the  behavior  of 
the  wave  is  determined  by  two  factors,  the  positive  impulse  y (which  is  conserved)  and 
y 

the  initial  slope It  can  be  shown  that,  to  a first  approximation,  the  "U"  wave  goes  over 

Tr 

to  an  N wave  whose  impulse,  y equals  the  impulse  of  the  first  positive  part  of  the  "U" 

V 

wave  Uu  of  Eq.  7-7)  and  for  which  - — is  equal  to  the  slope  of  the  "U"  wave  at  r = t,. 


Now 


dU 

dr 


V,. 


but  since  at  r = rr,  log 


T + T, 


r - r, 


7TT 

2t, 


log 


T + T(. 


2 rr, 


(7-17) 


, we  have 
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T(. 


0T 


IV  r 2 K 

z ‘ - - 2.5025574  -- 

wt,(t2  — r 2)  T> 


while  from  Eq.  (7-7b) 


(7- 1 8a) 


(7-1 8b) 


V,.t,  = .76373914  Vj, 

Solving  Eq.  (7- 18a)  and  Eq.  (7-18b),  we  get 

Vc  = 1.38249812  Vc  (7-19a) 

r,  = .55243413  t,  (7-19b) 

Hence,  the  equivalent  N wave  starts  out  with  an  amplitude  about  1.38  times  that  of  the  or.g.nal 

N wave  and  is  about  a little  more  than  half  as  long.  The  equivalent  N wave  is  shown  in  Fig  7- 

3.  To  this  approximation  the  peak  amplitude  is  given  by 


V 


peak 


.7637  VtT, 
.3996^-  + ( 


(7-20) 


and  the  period  of  the  positive  phase  by 


Asymptotically 


r„  = .8739 


U pe*  = -8739 


ir  = .8739 


(7-21) 

(7.22a) 

(7.22b) 


Thai  is.  both  (he  amplitude  and  length  are  .8739  of  what  they  would  have  been  without  the 
phase  shift,  i.e.  had  there  been  no  eaustte.  f-om  Eq.  (7-20)  is  compared  with  the  exact 

result  in  Table  I. 


C.  Improved  Analytical  Model  for  "U"  Wave 

The  disagreement  between  the  exact  result  and  the  equivalent  N wave  result  can  be 
explained  fairly  easily.  Consider  the  bottom  of  Fig.  7-2.  The  region  to  the  right  of  the  shock 
very  quickly  becomes  half  of  an  N wave  as  was  assumed  in  deriving  Eq.  (7-20).  The  deriva- 
tion, however,  ignores  the  area  to  the  left  of  the  shock  front.  This  area  is  not  negligible  until  f 
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--  -2 


Fig.  7.3  — The  asymptotic  form  of  Hie  equivalent  .V  wave  shown 
will  be  the  same  as  the  asymptotic  form  of  the  {/wave  for  large 
age 

is  quite  large.  Let  A,  be  the  area  to  the  left  of  the  shock  front  and  AR  be  the  area  to  the  right 
of  the  shock  front.  The  law  of  equal  areas  gives 

AR  + A,  = /t,  = (.7637)  yy  (7-23) 

If  we  assume  that  the  waveform  to  the  right  of  the  shock  is  an  N wave 

Ar-  y VU  + -3996  y (7-24) 

We  determine  AL  by  assuming  that  the  relatively  low  amplitude  waveform  to  the  left  of  the 
shock  is  undistorted.  The  retarded  time  of  the  shock,  is  given  by 


rs  = - VAZ  - .3996  y V,  - rr  (7-25) 

Since  UH  = U(t ,),  and  the  waveform  to  the  left  of  the  shock  is  assumed  to  be  undistorted, 
from  Eq.  (7-6) 


For  t J t,  » 1 we  can  use  the  approximation 


1 + tJt, 
1 - tJt, 


2 T, 
Ts  3 Tj 


(7-27) 


so  Eq.  (7-26)  becomes 


2Kt? 
3itt  , 


2 Vs* 


3rr \VA  f + ,3996y-  + r, 


Since  At  is  a correction  term,  we  let  VA  in  Eq.  (7-28)  be  given  by  Eq.  (7-20).  Substitut- 
ing the  resulting  expression  in  Eq.  (7-28)  and  using  Eqs.  (7-23)  and  (7-24),  we  obtain  for  the 
peak  value  of  the  Blokhintsev  variable 


.7637  y,.rt.  - - 

3 TT 


.7637—  ((  + .3996)  + .83356 

T,. 


I + 3996  y- 


(7-29) 


Equation  (7-29)  is  tabulated  in  Table  I.  It  is  seen  to  be  a much  better  approximation  when 
f > 1. 


D.  Solution  for  Linear  Case 


For  the  case  where  the  Gol'dberg  number  has  dropped  below  17.5,  the  solution  for  V 
takes  the  form  of  Eq.  (6-15)  rather  than  an  N wave.  Since  we  are  not  especially  interested  in 
the  exact  value  of  Fat  Z,.  but  rather  the  value  at  the  ground,  we  shall  use  the  same  approach 
for  the  linear  case  as  we  did  for  the  N wave.  That  is,  we  compute  the  linear  age  ( -E  right 
through  the  caustic,  as  we  did  in  Section  6,  and  find  the  waveform  by  applying  a w/2  phase 
shift  to  the  waveform  of  Eq.  (6-15).  The  phase  shift  does  not  affect  the  energy  or  spectral  con- 

£ 

tent  of  the  wave  so  V k and  — are  still  given  by  Eqs.  (6-21)  and  (6-22).  Since  will 


• X( 


■■■■II  I—  I I »l 


always  be  fairly  large  by  the  time  the  wave  gets  to  the  caustic,  we  compute  the  waveform  by 
applying  the  tt/2  phase  shift  to  the  asymptotic  form  for  V given  by  Eq.  (6-17),  The  result  is 


F(£,,  f) 


(— 1)*  A-! 
(2k  + D! 


(7-30) 


The  waveform  before  and  after  the  caustic  is  shown  for  £ / = 2 on  Fig.  7-4.  The  waveform  has 
two  small  positive  peaks  and  one  large  negative  one  given  by 


The  peak  pressure  at  the  ground  will  thus  be  negative  and  1.31  limes  larger  than  it  would  have 


been  had  there  been  no  caustic. 


I ig  7 4 — Typical  linear  wave  form  before 
and  after  rr/2  phase  shift  at  the  caustic 


8.  Results 


The  calculated  ray  paths  for  the  initially  upward  and  initially  downward  rays  (<f>  = 0°and 
<j>  = 180°,  i.e.  along  flight  path)  are  plotted  in  Fig.  8-1.  As  mentioned  in  Section  4,  the  rays  are 
straight  lines  below  Z = ZA  and  arcs  of  circles  [in  the  case  for  M = 2,  they  have  a radius  of 
2 (Z,  - Zt)]  between  Z = Z*  and  the  turning  point  at  Z = Z,.  The  formulas  for  computing 
the  range  for  the  various  segments  are  shown  on  the  figure.  The  initially  upward  ray  travels  a 
total  horizontal  distance  of  310  km.  The  initially  downward  ray  travels  about  20  km  farther. 
The  upward  traveling  wave  arrives  at  the  ground  in  22  minutes,  which  is  13  minutes  after  the 
Concorde  has  passed  overhead  (assuming  it  maintains  its  cruise  velocity).  That  is,  the  average 
horizontal  speed  of  the  sonic  boom  is  only  about  2/5  that  of  the  aircraft  at  cruise  conditions. 


sound  speed  (m/sec) 

0 300  600 


2?a.  »|  * 2^/3  (Zt-Zk) 4* — 

s/3  V 3 73 


Fig  8 I — Upward  and  downward  ray  paths  along  ground  track  lor  Concorde  flying  at  \l 
«*  2 using  our  assumed  sound  speed  profile  which  is  shown  on  the  right  of  the  figure 
Formulas  lor  computing  ray  distances  at  M - 2 are  shown  at  bottom  ol  the  figure 


The  initial  conditions  for  N wave  signature  are  found  using  the  formulas  of  Carlson.59 


The  formulas  derived  assume  that  the  asymptotic  form  of  the  wave  has  been  reached  and  the 
signature  can  be  characterized  by  the  bow  shock  overpressure  and  the  signature  duration. 
These  can  be  computed  in  terms  of  a shape  factor  (which  depends  on  the  aircraft  design),  the 
distance  from  the  aircraft,  and  the  length  of  the  aircraft. 
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Subscript  A refers  to  values  at  the  aircraft  altitude.  Values  of  AC,,  K,,  and  Kp  are  taken  from 
charts  in  Reference  39.  The  value  of  Af,  is  a function  of  the  lift  parameter  which  in  turn  is  a 
function  of  the  azimuthal  angle  (proportional  to  -cos  <b) . Thus,  the  lift  parameter  is  negative 
above  the  aircraft.  Carlson,  however,  does  not  compute  the  shape  factor  for  negative  lift 
parameters.  For  the  Concorde,  at  cruise  conditions,  the  shape  factor,  determined  by  an  extra- 
polation of  Carlson’s  figure  to  negative  lift  parameters,  would  have  a negative  value.  Since  we 
do  not  know  what  the  shape  factor  is  for  negative  lift  parameters,  we  have  assumed  a zero 
value  above  the  aircraft,  thus  neglecting  the  negative  lift  contribution  to  the  sonic  boom  in  the 
upward  direction.  This  will  result  in  an  overestimate  of  the  strength  of  the  upward  boom.  In 
the  downward  direction  we  can  accurately  obtain  the  shape  factor  as  a function  of  azimuthal 
angle.  We  approximate  Carlson's  figure  by  Ks  - Kn(A  — .6  cos  <A)  to  obtain  values  off  the 
flight  track. 

In  order  to  treat  the  boom  that  initially  starts  downward  and  is  then  reflected  from  the 
ocean,  we  first  calculate  the  boom  overpressure  and  duration  at  sea  level  using  Carlson's  for- 
mulas. In  order  to  compute  the  cylindrical  ray  tube  spreading,  we  start  the  problem  from  sea 
level  assuming  that  the  aircraft  is  at  the  mirror  image  point  below  sea  level  The  calculation  is 
then  carried  upward  to  the  top  and  then  back  down  again  For  calculations  of  the  upward 
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boom,  we  pick  a point  4 km  above  the  aircraft  to  start  the  problem.  For  the  off  track  calcula- 
tions, we  use  the  correction  formulas  of  Carlson  to  get  the  initial  signal  strength. 


The  results  for  <t>  = 0°,  and  <t>  = 180°,  for  the  initial  conditions  discussed  above,  are 
given  in  Figs.  8-2,  8-3,  8-4, and  8-5  for  the  initially  upward  wave  and  in  Figs.  8-6,  8-7  and  8-8 
for  the  initially  downward  wave  Figure  8-2  shows  Ap/p  as  a function  of  altitude  for  the  ini- 
tially upward  going  wave.  For  comparison,  the  results  of  a modified  "Gatwin's  Model",  dis- 
cussed in  Section  2,  are  also  shown.  For  this  comparison  the  Blokhintsev  invariant  was  held 
fixed,  while  the  sound  speed,  density,  and  ray  tube  area  were  varied  in  the  same  manner  as  for 
the  actual  calculations.  The  results  for  both  the  nonlinear-linear  (L)  and  nonlinear  only  (NL) 
solutions  are  orders  of  magnitude  smaller  than  Garwin's  model  predicts.  The  nonlinear  only 
solution  predicts  a maximum  value  of  Ap/p  of  .1,  which  corresponds  to  a Mach  number  of  .07. 
The  sharp  rise  in  A p/p  for  both  L and  NL  near  160  km  is  due  to  the  influence  of  the  caustic. 
From  Fig.  8-2  it  is  evident  that  once  the  linear  attenuation  takes  over,  the  Mach  number 
decreases  rather  quickly  In  Fig  8-3  the  Gol'dberg  number  is  plotted  as  a function  of  altitude. 
The  Gol  dberg  number  decreases  rapidly  with  altitude  in  the  vicinity  of  the  crossover,  as 
required  by  our  model.  Note  that  the  Gol'dberg  number  decreases  much  more  rapidly  after  the 
transition  value  when  linear  attenuation  is  included,  thus  justifying  our  "once  linear  always 
linear"  hypothesis  which  we  use  in  calculating  the  pressure  on  the  ground  for  cases  where  linear 
attenuation  becomes  important.  In  Figure  8-4  the  ratio  of  the  remaining  energy  to  the  initial 
energy  is  plotted  as  a function  of  altitude.  The  energy  decreases  rapidly  with  height  even  when 
linear  attenuation  is  not  included.  When  linear  attenuation  takes  over,  however,  the  decrease 
in  energy  is  much  more  marked.  Note  that  the  wave  has  lost  90%  ol  its  energy  by  the  time  it 
reaches  100  km  and  99%  of  its  energy  by  the  time  it  reaches  140  km.  Figure  8-5  shows  the 


fundamental  frequency 


2i„ 


which  is  seen  to  decrease  monotonically  with  altitude.  The  peak 


frequency,  is  2/3  of  the  value  shown 
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Fig.  8.2  — Relative  signal  strength  as  a function  of  altitude  for  an  intially  upward  go- 
ing ray  along  the  ground  track.  Nonlinear  only  INI.)  shown  by  solid  line.  Result 
with  linear  attentuation  included  (L)  shown  by  broken  lines  For  comparison,  the 
result  of  Garwin’s  model  is  shown  (no  attenuation  included). 


R GOL'DBERG  NUMBER 


Fig  8 3 — Gol'dberg  number  as  a fuction  of  altitude  This  number  is  a measure  of 
relative  importance  of  nonlinear  and  linear  effects  Solid  curve  (NL)  is  from  calcu- 
lation based  only  on  the  nonlinear  results 
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Fig.  8.4  — Ralio  of  energy  in  initially  upward  wave  to  initial  energy  as  a function  of  alti- 
tude. Again  NL  and  L refer  to  results  which  exclude  and  include  linear  attenuation 
respectively. 


FREQUENCY  (Hz) 


Fig.  8.5  — Fundamenial  frequency  of  initially  upward 
wave  (2/3  f^k)  as  a function  of  altitude 
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Figures  8-6  to  8-8  show  \p/p,  E/E„  and  fundamental  frequency  vs  altitude  for  the  ini- 
tially downward  wave.  The  initially  downward  wave  is  the  stronger  of  the  two  and  is  thus  the 
most  likely  signal  to  be  observed  at  the  ground.  This  is  compounded  by  the  fact  that,  due  to 
the  uncertainty  in  the  lift  factor,  we  have  probably  overestimated  the  initially  upward  boom. 
Due  to  the  higher  initial  amplitude  and  longer  length  of  the  initially  downward  wave,  the  transi- 
tion to  linear  attenuation  occurs  at  a higher  altitude  and  hence  linear  attenuation  is  a less 
important  factor.  Most  of  the  previous  observations  about  the  initially  upward  wave  hold  for 
the  initially  downward  wave.  Ninety  percent  of  the  energy  is  absorbed  below  100  km,  the  Mach 
number  never  exceeds  .1,  and  the  frequency  decreases  monotonically  with  altitude. 


Figure  8-9  shows  the  ground  pressure  for  the  initially  upward  and  initially  downward 
waves  as  a function  of  distance  from  the  ground  track.  The  solid  lines  show  the  combined 
linear-nonlinear  results,  the  dashed  lines  the  result  obtained  using  nonlinear  theory  alone.  The 
latter  can  be  considered  to  be  an  upper  bound  for  the  ground  pressure  The  dominant  signal  is 
from  the  initially  downward  wave.  The  highest  pressure  level  (about  .3  Pa)  will  occur  about 
300  km  from  the  ground  track.  The  pressure  measured  on  the  ground  track  is  a minimum  and 
will  be  about  .1  Pa  for  the  initially  downward  wave.  This  is  in  contrast  to  the  carpet  boom, 
where  the  pressure  is  a maximum  along  the  ground  track.  The  sharp  cutoffs,  indicated  by  the 
vertical  dashed  lines  at  380  and  425  km,  are  due  to  the  fact  that  for  values  of  1/cos  9„  < 1.13, 
the  waves  will  be  refracted  upwards  by  the  sound  speed  gradient  in  the  lower  troposphere 
before  reaching  the  ground. 

The  only  available  quantitative  data  directly  applicable  to  the  problem  of  Concorde  sonic 
boom  returns  from  the  thermosphere  has  been  taken  by  Lamonl-Doherty  Geological  Observa- 
tory/ They  have  bejn  monitoring  Concorde  sonic  booms  on  a regular  basis  using  a tripartite 
array  of  microbarographs  separated  by  about  1 km.  They  record  amplitudes  of  the  pressure  sig- 
nal with  either  a .1  to  1 Hz  or  1 to  10  Hz  bandpass.  We  compare  our  shock  strength  and  signal 
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Fig.  8.8  — Fundamental  frequency  of  initially  downward 
wave  (2/3  ) as  a function  of  altitude 
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Fig  8.9  — Relative  signal  strength  on  the  ground  from  a some  boom  returning  from 
the  thermosphere  as  a function  of  lateral  distance  Cutoff  is  determined  by  return- 
ing rays  which  are  refracted  upward  before  reaching  the  ground 
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period  with  their  peak  amplitude  and  frequency.  Lamont  reports  peak  pressure  of  1-5  Pa  and 
periods  in  excess  of  8 seconds  for  signals  identified  as  those  returning  from  the  thermosphere  at 
an  altitude  of  approximately  125  km  We  predict  an  amplitude  of  .3  Pa  and  a period  of  8.5 
seconds  for  a wave  returning  from  125  km.  Analysis  of  the  ray  tracing  indicates  that  the  origin 
of  this  wave  is  in  disagreement  with  the  analysis  given  in  their  report.  We  predict  the  origin  to 
be  further  away  and  at  bearing  1 1°  from  the  one  stated.  We  believe  the  origin  of  the  difference 
lies  in  the  difficulty  of  accurately  determining  bearing  at  such  low  frequencies  with  a 1 km 
array.  Lamont  lies  about  230  km  off  the  ground  track  we  believe  the  Concorde  to  be  on  while 
generating  the  signal  received.  From  Figure  8-9  we  see  for  this  distance  off  track,  the  received 
signal  would  be  as  large  as  will  be  measured  anywhere. 

Our  results  are  also  in  qualitative  agreement  with  the  measurements  done  at  Blacknest.40 
The  period  of  the  positive  phase  of  explosion  waves  which  return  from  the  thermosphere  are 
also  in  qualitative  agreement  with  the  results  shown  here,  which  indicates  that  our  theory  rea- 
sonably describes  shock  waves  which  propagate  to  the  thermosphere  and  return. 


9.  Strong  Shock,  Gravity,  Turbulence  and  MHD  Effects 


In  this  section  we  discuss  a number  cf  physical  effects  which  have  not  been  included  in 
our  model.  Our  objective  is  to  demonstrate  that  these  effects  will  either  have  no  effect  on  the 
solution  or  result  in  our  model  over  estimating  the  ground  pressure  levels 

A.  Strong  Shock  Effects 

The  intensification  of  upward  propagating  shock  waves  due  to  the  decreasing  density  is  a 
well  known  phenomenon.42  It  does  not,  however,  appear  to  be  an  important  consideration  for 
the  Concorde  problem.  Although  our  model  does  predict  substantial  increases  in  the  Mach 
number,  the  process  seems  to  be  self-limiting.  In  no  case  does  the  Mach  number  greatly 
exceed  0.1  which  is  a conservative  limit  for  the  validity  of  weak-shock  theory.4’  44  This  is  most 
likely  due  to  the  fact  that  the  wave  spends  proportionately  more  time  in  regions  of  lower  den- 
sity and  the  resulting  augmentation  of  the  aging  process  decreases  the  Mach  number  almost  as 
fast  as  the  decrease  in  density  increases  it.  The  smallness  of  the  Mach  number  assures  both  the 
validity  of  the  weak  shock  theory  and  that  the  wave  will  not  produce  any  serious  perturbations 
of  the  atmosphere. 

B.  Gravity  Effects 

When  an  element  of  fluid  is  displaced  upward  or  downward  from  its  equilibrium  position, 
the  gravitational  or  buoyancy  force  will  act  to  return  the  element  to  its  equilibrium  altitude  A 
fluid  element  thus  displaced  will  tend  to  oscillate  harmonically  with  an  angular  frequency  given 
by 
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This  frequency  is  known  as  the  Brunt-Vaisala  frequency.45  The  effect  of  this  additional  restor- 
ing force  is  to  introduce  dispersion  and  anisotropy  to  the  sound  speed  as  well  as  the  possibility 
of  a resonance  in  the  vicinity  of  the  Brunt-Vaisala  frequency.  None  of  these  effects  will  be 
important  as  long  as  the  frequencies  of  interest  are  large  in  comparison  with  the  Brunt-Vaisala 
frequency  or,  equivalently,  as  long  as  the  wavelength  is  much  less  than  the  scale  height. 


The  gravitational  field  will  have  its  most  pronounced  effect  for  vertical  propagation.  For 
this  case  the  phase  velocity  is  given  by4'’ 


v 


(9-2) 


At  160  km,  //,  = 26  km  while  A = 4.8  km.  Hence  v/c  = 1.0007.  At  17  km,  W,  = 6.9 
km  while  A = .03  km.  Hence,  v/c  = 1.0000004.  The  amplitude  effects  will  be  of  a similar 
size.  We  conclude  that  gravitational  effects  are  of  no  consequence. 


Turbulent  Effects 

In  the  above  discussions  we  have  assumed  that  the  atmosphere  was  steady  and  stratified 
horizontally.  The  real  atmosphere  is  known  to  be  unsteady  and  turbulent.  The  scale  of  the 
inhomogeneities  ranges  from  large  scale  weather  phenomena  such  as  fronts  and  storm  systems, 
to  internal  waves,  to  turbulent  fluctuations  which  may  have  a scale  length  down  to  fractions  of 
an  inch.  Most  of  the  work  that  has  been  done  on  the  effect  of  turbulent  fluctuations  on  sonic 
booms  has  looked  at  the  effect  in  the  region  below  3000  ft.  in  the  earth's  turbulent  boundary 
layer.’’1'  41  The  primary  effects  that  can  be  noticed  on  the  boom  signature  consist  of  random 
fluctuations  in  signal  strength  behind  the  shock  and  an  increase  by  several  orders  of  magnitude 
in  me  thickness  of  the  shock  front. 

Pierce  and  Maglieri’*  have  shown  that  atmospheric  turbulence  can  cause  variations  in 
waveform  (from  rounded  to  spiked  with  corresponding  peak  amplitude  variations  of  ± 8 dB). 
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These  variations  occur  both  in  time  (i.e.  as  measured  at  one  point  from  two  identical  aircraft 
passing  overhead  at  a 5 sec  interval)  and  in  space  (i.e.  as  measured  by  several  microphones 
along  the  flight  path  of  a single  aircraft).  This  variability  makes  the  problem  difficult  to  charac- 
terize. Work  by  Crow,48  George,49  Cole  and  Friedman,50  have  attempted  to  explain  the 
anomalous  rise  times  and  pressure  flucuations  in  terms  of  scattering  theories.  Pierce  attempts 
to  explain  the  results  in  terms  of  a wave  front  folding  mechanism.47  George  and  Plotkin51 
developed  a theory  where  the  nonlinear  steepening  just  balances  the  scattering  to  determine  the 
anomalous  rise  time  of  the  shocks.  All  of  these  mechanisms  tend  to  enhance  the  dissipation  of 
the  wave  or  leave  it  unchanged.  There  are  three  reasons  for  this: 

1.  Pierce  and  Maglieri38  have  shown  that  the  spectra  of  the  spiked  and  rounded 
waveforms  are  virtually  identical  for  frequencies  up  to  about  twice  ,/peak.  Since  long  range  pro- 
pagation is  determined  primarily  by  the  lower  frequencies,  one  would  expect  distant  arrivals  of 
both  waveforms  to  be  nearly  the  same. 

2.  The  Gol’dberg  number  for  the  spike  of  the  spiked  waveform  will  be  much  lower  than 
for  the  wave  as  a whole  because  its  duration  is  much  shorter.  One  would  therefore  expect  that 
in  the  thermosphere,  where  the  attenuation  is  high,  the  propagation  of  the  spike  would  be  dom- 
inated by  linear  attenuation  and  thus  the  spike  would  soon  disappear. 

3.  Turbulence  increases  the  effective  acoustic  attenuation.  Thus  the  Gol'dberg  number 
will  be  lower  and  linear  attenuation  will  take  over  sooner. 

Thus  a theory  that  ignores  the  effects  of  turbulence  on  the  signal  will  give  an  upper  bound 
on  the  signal  strength.  More  importantly,  the  highest  levels  of  turbulence  are  near  the  ground 
where  the  paths  we  are  considering  spend  very  little  time.  Most  of  the  attenuation  occurs  at 
the  higher  altitudes  in  the  thermosphere  where  turbulence  plays  no  role 
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In  addition  to  turbulence,  which  may  alter  the  local  signal  shape,  unsteady  inhomo- 
geneities  in  the  atmosphere  may  be  response  for  another  effect  which  has  been  observed  for 
signals  returned  from  the  stratosphere.  That  is,  most  of  these  signals,  rather  than  having  a sin- 
gle cycle  corresponding  to  an  N wave,  seem  to  have  several  cycles  which  may  correspond  to 
several  N waves.  Returns  from  the  thermosphere,  however,  seem  to  be  relatively  clean  having 
only  the  single  cycle.5  4I.  This  multiple  cycle  result  may  be  due  to  the  multipathing  of  signals 
which  are  bent  over  along  different  paths,  depending  on  the  unsteady  component  of  the  ray 
velocity.  This  effect  should  be  most  pronounced  near  where  the  ray  turns  over.  A slight 
change  in  wind  or  sound  speed  could  result  in  the  ray  traveling  a higher  or  lower  path.  The  fact 
that  this  mult, pie  signal  seems  to  come  primarily  from  the  stratospheric  and  tropospheric 
regions,  but  not  from  the  thermosphere,  may  be  a good  indication  that  the  thermosphere  is 
relatively  steady.  The  thermospheric  rays  pass  through  the  region  of  unsteady  flow  at  such  a 
steep  angle  that  they  are  virtually  unaffected  by  this  variation.  Multipathing  may  result  in  some 
enhanced  signal  if  two  paths  happen  to  come  together  with  their  phases  matched,  but  there  will 
also  be  regions  where  the  signal  would  be  diminished.  Increases  in  strength  of  the  signal  by 
more  than  a factor  of  two  seems  highly  unlikely  since  it  would  require  more  than  two  paths  to 
simultaneously  arrive  at  a point  in  phase.  The  major  effect  of  such  multipathing  seems  to  pro- 
duce signals  which  are  spread  out  in  time  and  involve  several  cycles. 

Another  effect  not  included  in  the  model  is  that  due  to  steep  wind  and  temperature  shear 
layers  which  have  a scale  length  of  the  same  order  as  the  boom  wavelength.  In  these  regions 
geometric  acoustic  theory  is  not  valid.  If  such  regions  exist,  they  will  produce  both  a reflected 
and  a transmitted  wave  from  the  shear  layer  (unless  the  velocity  becomes  targe  enough  to 
reflect  the  wave  completely).  If  such  shear  layers  are  weak,  they  will  have  only  a small  effect 
on  the  results,  but  in  any  case,  the  result  will  be  to  reflect  some  part  of  the  wave  energy  at  a 
lower  altitude.  Thus  our  model  will  provide  an  upper  bound  on  the  energy  which  reaches  the 
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thermosphere.  Any  energy  which  is  reflected  from  a lower  altitude  will  have  the  properties  of  a 
wave  which  would  have  been  refracted  at  that  altitude  (i.e.  will  have  a shorter  pulse  duration 
than  a wave  refracted  from  a higher  altitude).  Since  this  mechanism  will  tend  to  split  up  the 
energy  in  the  wave,  its  primary  effect  should  be  to  weaken  the  wave  We  are  not  aware  of  any 
evidence  which  suggests  that  any  strong  wave  is  reflected  to  the  ground  by  this  mechanism. 


Plasma  effects 


Since  the  atmosphere  at  160  km  is  partially  ionized,  the  possibility  of  exciting  MHD 
waves  which  may  propagate  at  a faster  rate  than  the  neutral  fluid  wave  exists.  If  we  look  at  the 
collision  frequency  for  scattering  of  charged  particles  by  the  neutral  fluid,  we  have 

IT  1/2 

— — sec1,  where  nt  is  the  neutral  density  and  a,,,,  is  the  cross  section,  approxi- 

m„ 

mately  equal  to  5 x 10'l5cm2.  A temperature  of  103  °K  and  a molecular  weight  of  about  20 
give  a collision  frequency  of  about  = 4.33  sec”1  or  about  4 collisions  per  second.  The  ratio 
of  ion  to  neutral  density  is  about  10  5 so  that  every  one  fourth  second,  10  5 of  the  energy  of 
the  wave  may  be  transferred  to  the  ions.  The  wave  spends  no  more  than  about  250  seconds  in 
the  region  above  140  km  (where  the  ion-neutral  ratio  is  >10”6).  Thus  about 
4 x 250  x 10 '5  = 10  2 of  the  neutral  energy  may  be  transferred  to  the  ions  as  the  wave  passes 
through  the  thermosphere.  Even  if  we  assume  that  this  energy  does  the  most  destructive  thing 
possible,  this  is  too  small  an  amount  of  energy  to  seriously  affect  our  results. 
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10.  Summary  and  Conclusions 


We  have  developed  a model  for  the  thermospheric  propagation  of  the  sonic  boom  from  a 
Concorde  aircraft  The  model,  which  is  for  the  most  part  analytic,  accounts  for  all  the  physical 
effects  which  we  believe  to  produce  a significant  effect  on  the  result.  From  the  results  obtained 
with  this  model,  we  reach  the  following  conclusions: 

1.  The  A pip  in  the  thermosphere  never  exceeds  .2.  This  implies  that  the  weak-shock 
model  is  valid  and  that  no  large  perturbations  of  the  thermosphere  (such  as  high  temperatures 
which  might  produce  strange  lights)  will  occur. 

2.  Most  of  the  energy  in  the  wave  is  lost  below  100  km.  There  is  thus  insufficient  energy 
remaining  in  the  wave  to  cause  winds  nearly  as  large  as  those  which  naturally  occur  in  the  ther- 
mosphere. 

3.  The  wave  will  reach  the  ground  at  a horizontal  distance  of  about  320  km  from  where  it 
was  launched. 

4 For  ray  paths  along  the  flight  track,  the  Concorde  will  pass  overhead  approximately  13 
minutes  before  the  first  arrival  of  the  boom  from  the  thermosphere  (assuming  that  the  Con- 
corde continues  at  its  cruise  speed  from  the  time  the  wave  is  launched).  The  average  horizon- 
tal velocity  of  the  sonic  boom  is  less  than  half  that  of  the  Concorde  at  cruise  conditions.  It  is 
thus  highly  unlikely  that  there  could  exist  reasonable  conditions  of  deceleration  under  which 
the  boom  would  precede  the  aircraft. 
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5.  The  sound  received  along  the  ground  track  is  a minimum  The  largest  signals  will  be 
observed  from  200  to  400  km  from  the  ground  track. 

6.  The  sound  received  will  in  all  cases  be  small  ( < .5  Pa)  and  of  very  low  frequency  K 
.1  Hz). 

We  conclude  that  thermospheric  returns  from  Concorde  are  of  sufficiently  low  amplitude 
and  frequency  that  it  is  unlikely  that  they  are  either  responsible  for  the  East  Coast  events  or 
likely  to  disturb  the  public.  In  addition  we  conclude  that  the  sonic  boom  from  the  Concorde 
does  not  have  sufficient  amplitude  or  energy  to  produce  a deleterious  effect  on  the  thermo- 
sphere. 
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